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1  EXECUTIVE  SUMMARY 

This  report  describes  the  research  performed  under  DTRA  grant  HDTRA 1-07- 1-00 16, 
“Thermal  History  using  Microparticle  Trap  Luminescence”  under  Principal  Investigators  Joseph 
Talghader  of  the  University  of  Minnesota  and  Eduardo  Yukihara  of  Oklahoma  State  University 
during  the  period  of  December  1 , 2006  -  May  31,2010. 

The  primary  objective  of  the  grant  was  to  investigate  a  thermal  sensing  technology  and 
associated  materials  that  extract  thermal  profiles  from  explosions  in  agent  defeat  tests  using 
micro-  and  nanoparticle  trap  luminescence.  The  number  of  occupied  traps  in  a  ceramic  varies 
dynamically  with  temperature.  Therefore,  if  the  traps  in  a  microparticle  are  filled  prior  to  an 
agent  defeat  test  using  deep  UV  light  or  ionizing  radiation,  then  the  trap  luminescence  afterwards 
is  a  probe  of  the  thermal  history  of  the  test. 

The  program  was  successful  in  all  of  its  major  components:  concept  demonstration, 
materials  research,  and  preliminary  explosive  testing.  The  laboratory  simulation  of  explosive 
heating  was  performed  using  micromachined  heaters  capable  of  reaching  hundreds  of  degrees  in 
milliseconds  or  less  and  cooling  back  to  room  temperature  with  arbitrary  profiles.  Among  the 
successful  concept  demonstrations  was  one  using  Mg2Si04:Tb,Co  particles  with  two 
thermoluminescent  peaks  heated  over  a  232  °C  to  313  °C  range  on  time  scales  of  less  than  200 
ms.  The  effect  of  maximum  temperature  during  excitation  on  the  intensity  ratio  of  the  two 
luminescent  peaks  was  compared  with  first -order  kinetics  theory  and  shown  to  match  within  an 
average  error  of  4.4%. 

In  materials,  lithium  fluoride  and  natural  topaz  were  identified  as  potential  temperature 
sensing  luminescent  materials.  A  new  luminescent  ceramic,  Mg2Si(>4:Tb,Co  was  also  found. 
This  material  has  two  traps  of  widely  different  energy  depths  within  the  band  gap  (with  a  weaker 
central  third  level).  TL  measurements  show  it  can  determine  peak  temperature  over  a  range  of 
roughly  150°C  to  400°C.  However,  it  appears  susceptible  to  optical  bleaching.  Oxyorthosilicate 
nanophosphors  were  also  investigated. 

Preliminary  explosive  testing  was  performed  in  collaboration  with  Jim  Lightstone  at 
Indian  Head  Naval  Surface  Warfare  Center.  Thermoluminescent  (TL)  A^OjiC  microparticle 
sensors  were  demonstrated  to  survive  detonation  of  pentaerythritol  tetranitrate  (PETN)  high 
explosives.  SEM  imaging  revealed  only  slight  shape  modification  during  an  explosive  event. 
Chemical  analyses  showed  small  increases  in  carbon  and  decreases  in  aluminum  indicating  A1 
combustion.  Finally,  the  TL  emission  peak  and  shape  remain  stable  to  within  1.9%  and  4.0% 
respectively,  less  than  experimental  error.  This  means  that  the  ability  of  the  particles  to  preserve 
thermal  history  information,  were  largely  unaffected  by  detonation. 
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2  INTRODUCTION 

When  a  cache  of  biological  weapons  is  destroyed,  every  region  containing  infectious 
particles  must  be  brought  to  a  minimum  thermal  profile  to  insure  sterilization.  Unfortunately,  the 
extreme  violence  of  an  agent  defeat  simulation  prohibits  existing  temperature  sensors  from 
functioning.  In  this  report,  the  investigation  of  a  thermal  sensing  technology  is  described  that 
extracts  thermal  profiles  from  the  luminescence  of  ceramic  microparticles,  particularly  oxides 
with  high  band  gaps.  The  number  of  occupied  traps  in  a  ceramic  varies  dynamically  with 
temperature.  Therefore,  if  the  traps  in  a  microparticle  are  filled  prior  to  an  agent  defeat  test  using 
deep  UV  light,  then  the  trap  luminescence  afterwards  is  a  probe  of  the  thermal  history  of  the  test. 

This  approach,  which  is  schematically  shown  in  Figure  1,  has  several  advantages: 

a)  The  particles  themselves  have  no  parts  that  can  be  damaged  by  the  agent  defeat  test; 
they  are  merely  implanted  quartz  or  similarly  tough  ceramics  that  can  easily 
withstand  the  harshest  environments  and  procedures. 

b)  The  materials  are  usually  very  cheap. 
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c)  The  particles  can  be  micromachined  to  take  on  just  about  any  shape  or  size,  and  thus 
will  have  essentially  identical  aerodynamics  and  thermal  experience  to  the  bioagents 
that  they  simulate. 

d)  Fluctuations  in  the  thermal  history  can  be  seen  from  the  differences  in  luminescence 
from  region  to  region  or  even  particle  to  particle. 

e)  The  luminescence  data  is  exponentially  weighted  towards  the  highest  temperatures 
experienced,  which  are  the  most  critical  in  destroying  infectious  agents. 

f)  Extreme  thermal  profiles,  such  as  those  seen  in  agent  defeat  tests,  can  be  simulated 
using  micromachined  heaters  that  increase  temperature  to  hundreds  of  degrees  and 
cool  in  less  than  a  millisecond. 

g)  The  luminescence  of  the  particles  is  largely  independent  of  particle  size  unlike 
methods  based  on  quantum  size  effects  or  phase  change.  This  means  that  precise 
control  of  size  is  not  critical  during  materials  fabrication. 


I.  Substrate  preparation 
1111  1111  l°n  implantation 
Quartz 

Silicon 


II.  Micromachining  to  make  particles 

Mask 


III.  Deep  UV  light  treatment 
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Figure  1  -  Measuring  the  thermal  history  of  an  agent  defeat  test  using  microparticle  luminescence.  The 
ceramic  particles  are  micromachined  to  match  the  size  and  shape  of  bacterial  or  viral  agents  and 
dispersed  in  a  burst  vessel.  After  the  test,  luminescence  from  the  microparticles  is  measured  to  extract 
the  temperature  experienced  in  any  region  of  the  system. 


3  METHODOLOGY 

3.1  Thermoluminescence  Testing  and  Apparatus 

The  simulation  of  explosive  heating  without  using  explosives  and  the  analysis  of  particles 
after  an  agent  defeat  test  both  required  a  new  method  of  luminescent  analysis.  Microheaters  were 
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ideal  for  these  as  they  were  easy  to  fabricate, 
had  extremely  short  time  constants,  and 
could  be  ramped  to  any  temperature  below 
their  degradation  points  with  the  application 
of  a  relatively  small  current. 

There  could  be  commercial  dual-use 
possibilities  for  microheater-driven  TL 
testing  as  well.  Many  radiation  dosimetry 
measurements  currently  use  badges  that 
contain  AbC^iC  particles  whose  trap  levels 
fill  based  on  the  radiation  dose  received. 
These  badges  are  then  periodically  sent  to  a 
central  laboratory  to  be  tested.  If  the  A^C^rC 
could  be  heated  in  pulses  with  microheaters 
with  an  integrated  avalanche  photodiode,  one 
might  be  able  to  obtain  real-time  analysis  of 
radiation  dose.  This  is  the  subject  of  a  recent 
patent  disclosure  of  the  UM/OSU 
collaboration. 

Figure  2a  is  a  scanning  electron 
microscope  (SEM)  image  of  a  row  of 
microheaters  from  early  in  the  program. 
They  were  fabricated  in  a  standard 
micromachining  process  using  a  polysilicon 
resistor  sandwiched  between  two  LPCVD 
nitride  layers.  Contact  to  the  poly  was  made 
using  an  RIE  etch  and  Cr/Au  metallization. 


Figure  3  -  Diagram  of  experimental  set-up.  The 
system  requires  very  little  complexity.  A 
photomultiplier  tube  (PMT)  is  placed  in  proximity 
to  the  microheater  and  connected  to  a  counter.  A 
current  source  drives  the  microheater. 


Figure  2  -  Images  of  early  microheaters,  a)  (top)  shows 
an  SEM  image  of  the  microheaters  after  fabrication,  b) 
(bottom)  shows  an  optical  micrograph  of  the  heaters  after 
AbOjiC  particles  have  been  dispersed  across  them. 

A  KOH  etch  was  used  for  release.  Figure  2b 
shows  an  optical  image  of  the  device  with 
A^C^C  particles  on  it. 

The  temperature  coefficients  of 
resistance  (TCR)  of  the  polysilicon  resistors 
were  determined  by  heating  the  chip  containing 
the  devices  using  a  thermoelectric  device  and 
measuring  the  resistance  changes  with 
temperature.  These  data  were  only  taken  over  a 
small  range  but  the  validity  of  the  TCR  data  to 
higher  temperatures  (~400°C)  was  strongly 
supported  by  the  linearity  of  the  1/R  vs.  I2  plot 
used  to  determine  the  microheater  thermal 
conductance. 

Figure  3  shows  a  diagram  of  an  early 
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experimental  setup.  After  Al203:C  particles  were  placed  on  the  heaters,  the  particles  were 
exposed  to  UV  radiation  to  fill  the  traps.  The  most  readily  accessible  source  was  a  deep-UV 
lamp  of  A,=205nm,  but  for  radioactive  testing,  a  p-decay  source  could  be  used.  After  irradiation, 
the  current  to  the  microheaters  was  gradually  increased.  We  estimated  a  heating  rate  of 
2°C/second.  As  the  particles  were  heated  they  emit  light  whose  intensity  is  related  to  the 
radiation  dose  they  have  received.  The  peak  of  the  TL  is  indicative  of  the  depth  of  the  traps 
within  the  bandgap. 


In  the  second  and  third  years,  we  completely 
automated  our  basic  thermoluminescence  (TL) 
measurement  system.  A  particle  of  TL  material  can 
now  be  placed  on  a  microheater,  and  once  the 
computer  control  is  activated,  the  traps  can  be  filled 
by  UV  light,  an  explosive  heating  cycle  can  be 
performed,  a  TL  ramp  can  be  performed,  and  then 
the  cycle  can  be  repeated  as  often  as  necessary  to 
obtain  the  full  relationship  between  thermal  history 
and  luminescence.  Figures  4,  5,  6,  and  7  below 
show  the  system  details. 

At  the  core  of  our  thermoluminescence  curve 
collection  setup  is  a  microheater  (Figure  4) 
consisting  of  a  150nm-thick  platinum  resistor  Figure  4  -  A  microstructure  used  for  heating 

sputtered  on  lnm  adhesion  layer  of  chrome,  in  turn  the  luminescent  microparticles  at  rates 

Trvn  t  j  .  ...  simulating  agent  defeat  explosions, 

on  a  200nm  LPCVD  nitride;  silicon  wet-etching 

makes  this  a  platform  suspended  above  a  ~300pm- 

deep  pit.  The  linear  temperature  vs.  resistance  relationship  of  the  resistor  is  found  by  measuring 
its  resistance  while  heating  it  on  a  hotplate.  The  microheater  is  driven  by  a  Keithley  2410 
sourcemeter  for  linear  heating  ramps  or  a  HP  8114A  pulse  generator  for  explosive  heating. 
Particle  luminescence  is  measured  by  a  Electron  Tubes  Ltd.  P25PC  photomultiplier  tube  and  BK 

Precision  1823A  frequency 
counter;  the  radiation  for  filling 
sample  traps  is  provided  by  a 
deuterium  lamp  putting  out  0.8 
mW/cm2  of  180-400nm  UV  light. 
A  motorized  aluminum  shutter 
provides  protection  for  the  PMT 
during  irradiation.  (Figures  5-7). 

The  script  for  an  explosively- 
heated  thermoluminescence 
experiment  is  to  first  fill  the  traps 
via  UV  irradiation,  then 
explosively  heat  the  sample  via 
voltage  pulse  to  the  microheater. 
A  linear  heating  ramp  then 
follows  to  induce  and  measure 
thermoluminescence.  For 


PC 


FREQUENCY  COUNTER 


MONOCHROMATOR 


PHOTOMULTIPLIER 


SHUTTER 


PULSE  GENERATOR 

1 

Supplying  voltage  pulses 

PARTICLE(S) 

W 

MICROHEATER 


Figure  5  -  Conceptual  diagram  of  experimental  apparatus. 
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normalization  and 

consistency  checking 

purposes,  a  non-pulsed  TL 
reference  curve  is 
sometimes  collected  before 
the  traps  are  refilled  for  the 
main  pulsed  experiment. 

A  number  of  custom 
Python  programs  coordinate 
and  control  most  of  the 
thermoluminescence 
measurement  process.  This 
allows  the  experiment  to 
run,  automated  and 
remotely-controlled,  inside  a 
blacked-out  darkroom 

admitting  an  average 

measured  background  light  signal  of  around  600  counts  per  second  (including  the 
photomultiplier's  quoted  200  CPS  dark  count  rate).  The  recent  addition  of  feedback  sensing  to 
many  experiment  parts  provides  enhanced  stability  and  confidence  of  precision.  Resulting  data  is 
processed  mainly  using  MATLAB. 


1.  XeBWey  2410  souroenrwrtw  drives  mctcDeators  Tor  nneat  ramp 
0  Acjlent  HkxiM  DC  *upply  axnpoii*  ctxrfrotocf  pnwrx  lot  PMT 

3  SK  Precision  1823A  frequency  counter  /roots  PMI  counts 

4  HoathM  TP-2713  DC  suotfy  powers  stopper  control  cttcutr 

5.  HP  3H4A  pu&c  gorcicfcx  Tor  microtKOtor  oxpiosrw  heating 


6  W  source*  s-ves  "8d-6GQ  nm  light  for  sampte  modiatton 

7.  If D  dnitooom  bulb  /cv  txmdhng  kQM-sansttiYH  samfitos 

8.  Mtefoprobe  po*boner*  cornea  to  mrcrohoator  contacts 

9.  pc  controls  f-6  aeovs  for  automateryremoto  operation 


Figure  6  -  Photograph  of  the  control  electronics  of 
thermoluminescence  excitation,  and  testing  system. 


the 


10.  Custom  sample  chamber  holds  UV  source,  PMT,  and  shutter  assembly  in  place 

11.  Stepper  motor  moves  aluminum  shutter  (inside  chamber)  with  position  sensing 

12.  Stepper  motor  control  circuit  translates  PC  commands  to  shutter  movement 

13.  PC  power  supply  provides  high  currents  required  by  stepper  motor 

14.  Electron  Tubes  P25PC  PMT  high-sensitivity  light  detection  for  thermoluminescence 


Figure  7  -  Shuttering  system  and  sample  chamber  for  microheaters. 

Obviously,  one  cannot  perform  actual  agent-defeat  testing  in  a  university  setting.  In  order 
to  eliminate  the  need  to  explode  particles,  we  perform  much  of  our  thermoluminescent  testing  on 


HDTRA1  -07-1  -001 6 

University  of  Minnesota  (Twin  Cities) 


micromachined  heaters.  These 
devices  are  fabricated  using  VLSI 
techniques  and  are  essentially 
dielectric  plates  with  integrated 
metal  or  semiconductor  resistors. 
The  time  constants  of  these  heaters 
vary  from  several  tens  of 
microseconds  to  hundreds  of 
milliseconds  depending  on  their  size 
and  thermal  isolation.  Figure  8 
shows  the  Temperature  versus  Input 
Power  for  a  heater  than  is  directly 
connected  to  its  substrate.  The  time 
constant  for  this  device  is  about 
1  OOps,  which  is  even  faster  than  one 
expects  in  the  periphery  of  an  agent 
defeat  event.  Figure  9  shows  an 
SEM  image  of  a  device  of  this  type. 
(Of  course,  in  the  center  of  an  agent 
defeat  event,  the  heat  transfer  time 


will  be  even  taster  than  this  because  the  transfer  is 
mediated  by  radiation;  however,  the  peak 
temperatures  in  a  radiation-limited  environment  are 
so  high  that  no  infectious  agents  would  survive 
anyway.) 

Figure  10  shows  the  Temperature  versus 
Power  curves  for  a  microheater  whose  substrate  has 
been  etched  out  from  underneath  it.  This  type  of 
microheater  can  be  heated  to  hundreds  of  degrees 
assuming  the  dielectrics  and  metals/semiconductors 
have  been  chosen  properly.  The  time  constant  of 
such  microheaters  is  generally  in  the  tens  of 
milliseconds,  which  is  sufficient  for  most  explosive 
heating  simulations  and  is  certainly  capable  of 
distinguishing  between  different  levels  of  a  slow 
cool-down. 

Figure  11  shows  a  potential  pitfall  of 
microheaters  if  they  are  not  designed  properly: 
nonuniform  heating.  This  problem  was  seen  in  the 
first  year  microheaters  but  was  subsequently  remedied  by  insuring  that  the  heating  resistors  ran 
through  every  portion  of  the  heaters.  Images  of  such  heaters  are  shown  in  Figure  9  and  in  the 
results  sections.  Another  potential  pitfall  is  ensuring  that  the  particle  thermal  mass  is  not  larger 
than  that  of  the  microheater.  Otherwise,  the  particle  temperature  will  not  follow  the  heater 
temperature. 


Figure  9  -  High  speed  microheater.  This 
heater  has  a  short  time  constant  because  it 
is  directly  connected  to  the  substrate; 
however,  it  cannot  reach  the  temperatures 
of  released  heaters  without  damage. 


Figure  8  -  Timing  diagram  of  a  microheater  similar  to  that 
shown  Figure  6.  Note  that  the  time  to  reach  275°C  is  less 
than  200ps. 


Normalized  Intensity 
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Membrane  Temperature  vs.  Power  Input 
Lot  1,  Wafer  11,  Device  4A-2 


Figure  10  -  Temperature  vs.  Power  for  a  released  microheater. 


Figure  11  -  Repeatability  and  variability  in  early  microheater  TL.  On  the  left,  Figure  (a)  shows  that  a 
microparticle  on  a  heater  that  is  excited  by  deep  UV  light  repeatedly  and  tested  repeatedly  in  the  same 
position.  It  has  extremely  stable  TL  characteristics.  On  the  right,  however,  particles  are  tested  on 
different  regions  of  a  microheater.  Particles  located  near  the  heating  resistors  on  the  microheaters  are 
apparently  exposed  to  higher  temperatures  than  particles  located  away  from  them.  Based  on  these 
results,  we  redesigned  our  heaters  to  have  resistors  that  covered  eveiy  region  on  the  heater  for 
temperature  uniformity  so  that  all  particles  gave  the  same  TL  regardless  of  location. 
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3.2  Samples 

This  project  included  the  characterization  of  existing  materials  (commercially  available  or 
obtained  through  collaborators),  and  experimental  characterization  of  materials  synthesized  in 
our  laboratory  at  Oklahoma  State  University.  Table  1  lists  the  materials  investigated  in  this 
project,  the  synthesis  process  used  to  produce  them,  and  the  source  of  the  material  (if  not 
synthesized  in  the  laboratory). 


Table  1,  List  of  materials  characterized  for  potential  use  as  temperature  sensors. _ 

Material _ Synthesis  process  Source _ 

Dosimetric  materials 

A1203:C  Crystal  growth  Landauerlnc. 

LiF:Mg,Ti  Crystal  growth  Thermo  Fisher  Scientific  Inc. 


Artificial  materials 

LiA102 

Sintering 

Mg2Si04 

Sintering 

CaSi03 

Sintering 

Y2Si05:Ce  (YSO) 

Crystal  Growth 

Gd2Si05:Ce  (GSO) 

Crystal  Growth 

Lu2Si05:Ce  (LSO) 

Crystal  Growth 

Y2Si05:Ce  (YSO) 

SCS 

Gd2Si05:Ce  (GSO) 

SCS 

Lu2Si05:Ce  (LSO) 

SCS 

A1203 

SCS 

Li2B407 

SCS 

MgO 

SCS 

La203 

SCS 

LaAI03 

SCS 

LaMg03 

SCS 

MgAI204 

SCS 

LaMg  Al  1 1 0 1 9 

SCS 

Ce02 

SCS 

Institute  of  Nuclear  Sciences,  Belgrade,  Serbia 
Institute  of  Nuclear  Sciences,  Belgrade,  Serbia 
Institute  of  Nuclear  Sciences,  Belgrade,  Serbia 
Los  Alamos  National  Laboratory 
Hitachi  Chemical  Co.,  Ltd. 

Los  Alamos  Natinal  Laboratory 
Los  Alamos  National  Laboratory 
Los  Alamos  National  Laboratory 
Los  Alamos  National  Laboratory 


Natural  materials 

Al2Si04(F,0H)2  (topaz)  Natural 

Al2SiQ5_ Natural 


3.3  Material  synthesis 

Samples  were  synthesized  at  Oklahoma  State  University  using  the  Solution  Combustion 
Synthesis  (SCS)  technique  (Chick  et  al.,  1990;  Kingsley  and  Patil,  1988;  Kingsley  et  ah,  1990; 
Yen  and  Weber,  2004).  The  reagents  were  dissolved  in  water,  dried  in  a  hot  plate  at  200°C  for 
~T-2  h.  The  temperature  of  the  hot  plate  was  then  increased  to  500°C  until  combustion  occurred 
(~7  min).  The  remaining  powder  was  collected,  crushed  in  a  agate  mortar,  and  annealed 
(typically  at  900°C  for  2  h)  to  remove  organic  residuals  and  promote  growth  of  the  crystallites. 
Doped  materials  were  investigated  by  adding  the  dopant  (e.g.  Ce(N03)3-6H20)  in  the  initial 
solution  in  the  desired  amount.  Figure  12  illustrates  the  stages  of  the  SCS  process  and  Table  2 
lists  the  materials  synthesized,  typical  quantities  used  in  the  synthesis,  and  dopants  investigated. 


HDTRA1  -07-1  -001 6 

University  of  Minnesota  (Twin  Cities) 


Figure  12.  Stages  of  the  Solution  Combustion  Synthesis  process.  The  first  photo  shows  the  material  after  the 
1  h  drying,  after  which  the  hot  plate  temperature  was  increased  to  500°C  and  the  other  photos  were  obtained. 
Combustion  is  seen  -6  min  after  raising  the  temperature  of  the  hot  plate. 


Table  2.  Materials  synthesized  by  SCS  and  dopants  investigated. 


Material 

Reagents 

Dopants 

A1203 

5  g  A1(N03)3-9H20, 2.0  g  (NH2)2CO  (urea); 

18.7  g  A1(N03)3-9H20,  6.25  g  C2H502N  (glycine) 

Ce,  Er,  Li 

Li2B407 

10.2  g  NH4NO3,  2.47  g  H3B03,  3.69  g  Li2C03,  7.65  g  (NH2)2CO 
(urea) 

Ce,  Eu 

MgO 

12.8  g  Mg(N03)2-6H20,  5  g  (NH2)2CO  (urea); 

12.8  g  Mg(N03)2-6H20, 4.166  g  C2H502N  (glycine); 

12.8  g  Mg(N03)2-6H20,  2.086  ml  C2H502N  (ethylenediamine); 

Ce,  Er,  Li,  Tm,  Yb,  Tb,  Eu 

La203 

21.6  g  La(N03)2-6H20,  7.5  g  (NH2)2CO  (urea) 

Ce,  Li 

LaA103 

LaMg03 

1 1 .5  g  La(N03)2-6H20,  10  g  A1(N03)3-9H20,  4  g  (NH2),CO  (urea) 

Ce 

MgAl204 

5g  A1(N03)3-9H20,  1.71  g  Mg(N03)2-6H20,  2.675  g  (NH2)2CO 
(urea); 

Ce,  Er 

LaMgAlnOi9 

2.1  g  A1(N03)3-9H20,  1.24  g  Mg(N03)2-6H20,  20  g 

A1(N03)3-9H20,  9.2  g  (NH2)2CO  (urea); 

Ce 

Ce02 

10.8  g  Ce(N03)2-6H20, 3.75  g  (NH2)2CO  (urea) 

Eu,  Li 
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3.4  Equipment 

3.4.1  Radioluminescence  (RL)  measurements 

RL  measurements  were  performed  using  a  Magnum  X-Ray  tube  (Moxtek,  Orem,  USA) 
operating  under  40  kVp  with  W  filament  and  Ag  target.  All  measurements  were  made  using 
100  pA  tube  current,  which  results  in  a  dose  rate  of  approximately  1 50  mGy/s.  The  RL  was 
detected  using  a  USB-2000  miniature  fiber  optic  spectrometer  (Ocean  Optics),  equipped  with  a 
ILX51 1  linear  silicon  CCD  array  (Sony)  and  a  fixed  polymer  rule  grating  allowing  a  700  nm 
dispersion  range  with  maximum  efficiency  at  500nm.  The  luminescence  is  guided  from  the 
sample  toward  the  spectrometer  using  an  optical  fiber  (1  mm  core  diameter,  transmission 
between  200  -  1 1 00  nm)  and  a  f/2  fused  silica  lens  to  couple  the  sample  luminescence  to  the 
optical  fiber.  The  spectrometer  has  an  order-sorting  detection  filter  to  eliminate  the  second  order 
interference  light  and  - 7.2  nm  resolution. 

3.4.2  Thermoluminescence  (TL) 

TL  measurements  were  performed  using  a  Riso  TL/OSL-DA-15  reader.  Irradiations  were  carried 
out  using  the  integrated  ^°Sr/90Y  beta  source  delivering  a  dose  rate  of -100  mGy/s.  The  TL  was 
detected  using  an  Electron  Tubes  9235QB  photomultiplier  tube  with  optical  filters  in  front  of  it 
to  select  specific  emission  bands.  Typical  filters  used  are  Schott  BG-39  for  visible  detection 
(6  mm  total  thickness,  Schott  AG,  Mainz,  Germany)  and  Hoya  U-340  filters  for  UV  detection 
(7.5  mm  total  thickness,  Hoya  Corporation,  Tokyo,  Japan  ).  Nitrogen  gas  was  purged  in  the 
sample  chamber  for  the  TL  measurements  to  avoid  spurious  luminescence. 


3.4.3  X-ray  diffraction 

X-ray  diffraction  data  was  obtained  using  a  Phillips  Analytical  X-ray  PW3020  with  Cu  Ka 
radiation  by  scanning  the  20 in  0.02  degrees  step  size  and  0.5  s  scan  step  time. 


4  RESULTS 

4.1  Modeling  the  Effect  of  Thermal  History  on  Idealized  Materials 

In  this  section,  we  briefly  examine  the  ability  of  ideal  first  order  materials  to  distinguish  time  and 
temperature  in  a  rapid  thermal  event.  The  next  section  will  discuss  real  luminescent  materials 
and  section  4.3  will  cover  actual  experimental  demonstration  of  temperature  sensing  in  rapid 
thermal  events. 

The  ultimate  goal  of  using  microparticle  thermoluminescence  for  thermal  history  sensing 
is  to  be  able  to  extract  the  entire  history,  not  just  maximum  temperature.  Further  one  wishes  to  be 
able  to  distinguish  extremely  rapid  changes.  The  first  step  towards  this  ultimate  goal  is  to 
independently  be  able  to  isolate  maximum  temperature  and  cooling  time.  (Note  that  heating  rate 
is  extremely  fast  in  explosions  and  therefore  plays  less  of  a  role  in  the  total  energy  delivered  to 
an  object  or  biological  particle  than  Tmax  and  cooling  rate.)  The  following  discussion  shows  that 
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a  variety  of  Tmax  and  cooling  combinations  can  give  identical  results  in  a  TL  curve.  For 
example,  consider  that  the  periphery  of  an  explosion  has  a  maximum  temperature  of  800K  and 
cools  over  10  seconds.  This  thermal  history  would  cause  a  shallow  trap  in  a  TL  material  to  empty 
but  a  deeper  trap  may  still  be  populated  this  leads  to  a  specific  population  ratio.  If  this  is  the  only 
pair  of  traps  used  then  an  explosion  with  (Tmax  =  800K,  cooling  =  10  seconds)  may  give  the  same 
TL  results  as  an  explosion  with  (Tmax  =  81  OK,  cooling  =  5  seconds).  However,  if  one  uses  three 
or  more  traps,  one  can  specify  both  Tmax  and  cooling  time  by  looking  at  the  data  from  two  or 
more  different  trap  population  ratios. 

A  single  first-order  trap--i.e.  one  combination  of  s  and  E— which  undergoes  an  explosive 
heating  event  of  a  given  maximum  temperature  and  duration  with  initial  electron  population  no 
may  afterward  host  a  certain  population  of  still-trapped  electrons.  Under  first-order  kinetics 
theory,  for  each  possible  value  of  one  thermal  event  parameter  (either  maximum  temperature  or 
cooling  rate),  there  is  a  single  value  of  the  other  which  will  result  in  exactly  the  same  surviving 
population.  If  the  difference  between  the  population  number  resulting  from  an  “actual”,  pre¬ 
determined  explosive  heating  event  and  those  from  another  set  of  explosion  parameters  are 
plotted  on  the  plane  of  explosive  heating  and  duration  variables,  the  result  is  a  one-to-one  curve 
along  which  the  match  is  exact,  indicating  that  the  explosive  heating  event  could  have  possessed 
any  pair  of  parameters  lying  along  this  “correctness  curve”.  Given  one  of  the  explosive  heating 
event’s  two  distinguishing  values,  the  other  can  be  found. 


The  correctness  curve,  shown  in  3D  and  contour  views  in  Figures  13  and  14,  were 
calculated  from  two  traps  subjected  to  a  simulated  explosive  heating  event  reaching  to  300  °C 
lasting  for  5 

10CU 


and 

seconds.  A  similar  but 
not  identical  profile  is 
shown  in  Figure  15. 
This  set  of  parameters 
can  indeed  be  found  on 
the  line  along  which  the 
correctness  is  100%,  i.e. 
where  the  actual 

remaining  population  of 
the  traps— the  only 
information  about  the 
explosion  that  the 
producing  simulation  is 
given— matches  that 
predicted  by  theory.  For 
visualization  purposes, 
we  have  found  it  easiest 
to  calculate  the  percent 
error  between  the 
“actual”  explosive 
heating  event  parameters 
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Figure  13  -  Correctness  curve  that  represents  how  identical  TL  signatures  may  be 
obtained  from  different  combinations  of  maximum  temperature  and  cooling.  The 
“correctness”  is  defined  by  how  well  the  population  ratio  of  two  peaks  (for  example 
the  two  MgSi04:TbCo  peaks  in  Figure  4)  matches  that  of  one  specific  Tmax/cooling 
time  combination.  Note  that  many  combinations  with  lead  to  identical  TL  results  if 
only  one  particle  or  pair  of  peaks  is  used. 
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and  the  candidates  being 
considered,  cap  large  error 
values  at  100%,  and  then  to 
subtract  these  from  unity  to 
yield  what  we  term  the 
“percent  correctness”. 

To  provide  a 
reference 

thermoluminescence  level 
with  which  to  counteract  the 
effects  of  varying  overall 
luminescence  from 

observation  to  observation, 
another  trap  may  be  added 
and  the  ratio  between  the  two 
glow  curve  peaks’  heights 
used  instead  of  the  absolute 
intensity.  Therefore,  two 
traps  of  suitably  differing 
parameters  will  in  theory 
enable  the  experimental 
recovery  of  a  single 
explosive  heating  parameter— 
either  maximum  temperature 
or  duration— if  the 

other  parameter  is 
already  known. 


Changing  the 
trap  parameters  causes 
the  shape  of  the 
correctness  curve  to 
change.  In  order  to 
reconstruct  the  thermal 
event  without  first 
knowing  either 

parameter,  we  can 
therefore  introduce 
additional  traps  to 
undergo  the  same 
heating;  if  the  trap 
characteristics  are 
different  enough,  the 
resulting  lines 

determining  the 

explosive  heating 


3  4  5  6  7 

Cooling  time  (sec) 

Figure  14  -  Contour  plot  of  “correctness  curve”  that  represents  how  identical 
TL  signatures  may  be  obtained  from  different  combinations  of  maximum 
temperature  and  cooling.  The  “correctness”  is  defined  by  how  well  the 
population  ratio  of  two  peaks  (for  example  the  two  MgSi04:TbCo  peaks  in 
Figure  4)  matches  that  of  one  specific  Tmax/cooling  time  combination.  Note 
that  many  combinations  with  lead  to  identical  TL  results  if  only  one  particle 
or  pair  of  peaks  is  used. 


Figure  15  -  Heating  profile  to  simulate  an  explosive  event.  A  rapid 
linear  heating  ramp  bring  the  temperature  to  400°C  in  a  few 
milliseconds.  Afterwards  a  slow  cooling  over  several  seconds  brings 
the  temperature  back  to  ambient. 
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parameters  that  each  trap’s  (or  pair  of  traps’)  results  deem  possible  intersect,  giving  us  a  single 
set  of  candidates  to  describe  the  explosive  heating. 


The  overlaid 
correctness  curves  of  three 
pairs  of  traps  subjected  to 
the  same  thermal  heating 
event  reaching  to  300°C 
and  lasting  for  5  seconds. 
Subsequent  calculation 
using  only  knowledge  of 
the  trap  parameters  and  the 
pre-  and  post-explosion 
electron  populations,  as 
can  be  gathered  by 
experiment,  reveals  that 
the  intersection  of  the 
three  correctness  curves 
indicates  a  single  possible 
thermal  event  description: 
a  maximum  temperature 
of  300°C  and  a  duration  of 
5  seconds.  By  using  two 
correctness  curves  and  a 
single  shared  reference 
trap,  the  same  result  can 
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Figure  16  -  Diagram  showing  two  contour  plots  from  two  pairs  of  traps  with 
similar  characteristics.  Note  that  even  these  relatively  similar  trap  pairs  give 
independent  curves  with  a  well-defined  intersection  that  can  identify  the  true 
maximum  temperature  and  cooling  time  of  an  agent-defeat  event. 


be  achieved  with  just  three  traps,  although  some  tuning  of  their  parameters  is  necessary. 


max 

(°C) 


ni/n2: 

n2/n3: 

=0.785 

=1.496 

J  1  VJ 

300  ni/n2= 

0.867 

k. 

"actual' 

^  c 

J 

scenario 

nj/n2=C 

.840 

300  n2/n3= 

9QH 

1.262  C 

)  k 

ni/n2 

n2/n3 

r 

7  ^ 

=0.853 

=1.295 

J  n2/n3=l 

.329 

l 

1 

K 

n1/n2: 

n2/nr 

V. 

J 

:0.902 

=1.177 

e 

\ 

Particles  used: 

1:  E=1.277  eV,  s=lxl012 
2:  E=1.30  eV,  s=lxl012 
3:  E=1.30  eV,  s=2x!012 


^ cooling  (sec) 

Figure  17  -  This  graph  shows  the  ratios  of  three  traps  with  different  characteristics  (trap  depth 
and  phonon  interaction  frequency,  shown  at  right)  under  the  influence  of  different  heating 
conditions.  Note  that  one  can  distinguish  both  time  and  temperature  under  these  conditions.  The 
traps  are  *not*  the  same  as  those  simulated  in  Figure  16. 
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4.2  Modeling  the  Effect  of  Thermal  History  on  Luminescence  Materials 

The  behavior  of  real  luminescence  materials  with  thermal  history  or  pre-heating  temperature  is 
dependent  not  only  on  the  maximum  temperature,  but  on  the  time-temperature  profile.  To 
determine  the  thermal  history  during  explosive  events,  it  is  necessary  to  develop  an 
understanding  on  how  much  the  heating  and  cooling  rates  affect  the  population  of  trapped 
charges. 

To  develop  this  understanding,  we  simulated  a  material  with  three  trapping  levels  at  different 
depths  below  the  conduction  band,  which  after  irradiation  produces  the  hypothetical  TL  curve 
shown  in  Figure  18.  Each  one  of  these  TL  peaks  provides  information  on  a  specific  temperature 
range. 

The  effect  of  pre-heating  on  the  TL  curve  is  shown  in  Figure  19a.  As  the  maximum  pre-heating 
temperature  increases,  the  intensity  of  the  TL  peaks  decrease  sequentially  because  the 
corresponding  trapped  charges  are  emptied.  The  intensity  of  the  TL  peaks  is  plotted  against  the 
maximum  temperature  in  Figure  19b.  The  range  in  which  the  peaks  decrease  is  the  useful 
temperature  range  that  can  be  “probed”  by  each  TL  trap. 


Step  Annealing  (Heating  Rate  IC/s),  Tph=43 


Figure  18.  TL  curve  of  a  hypothetical  material  containing  three  trapping  levels  at  depths  1.19  eV,  1.33  eV, 
and  1.53  eV.  The  frequency  factor  was  assumed  to  be  1012  s'1  and  the  heating  rate  was  l°C/s. 
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Figure  19.  Effect  of  maximum  pre-heating  temperature  on:  (a)  the  TL  curve  of  a  hypothetical  material  with 
three  trapping  levels;  (b)  the  maximum  intensity  of  the  TL  peaks. 

However,  the  “useful”  temperature  range  of  each  TL  trap  depends  in  general  on  the  heating  rate. 
Therefore,  we  simulated  how  the  traps  are  emptied  by  pre-heatings  applied  at  various  heating 
rates,  from  0.1  K/s  to  500  K/s.  The  trap  population  is  plotted  against  the  maximum  pre-heat 
temperature  in  Figure  20. 

As  the  heating  rate  of  the  thermal  history  increases,  the  detrapping  profile  shifts  to  higher 
temperatures  and  becomes  less  dependent  on  heating  rate.  Therefore,  in  the  case  of  a  fast  rise  in 
temperature  characteristic  of  an  explosion,  the  population  of  trapped  electrons  will  depend  only 
on  the  maximum  temperature,  and  not  on  the  heating  rate. 


z  *et<? 


Figure  20.  Trap  population  versus  maximum  pre-heat  temperature  for  heating  rates  from  0.1  K/s  up  to 
500  K/s. 

The  simulations  carried  out  so  far  indicated  that  the  possibility  of  determining  the  maximum 
thermal  history  temperature  is  more  feasible  for  fast  heating  rates  (explosion)  than  for  slow  ones, 
but  one  has  to  keep  in  mind  the  “useful”  temperature  range  of  each  trapping  center.  However,  the 
simulations  assume  an  instantaneous  cooling  to  room  temperature.  Influence  of  the  cooling  rate 
will  be  discussed  in  the  next  section. 
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4.2.1  Test  of  TL  models 

The  investigations  involved  in  the  development  of  temperature  sensing  using  TL  materials 
require  modeling  of  the  TL  process.  These  models  can  be  useful  to  understand  the  influence  of 
extremely  fast  heating  rates  or  to  predict  the  effect  of  different  temperature  profiles  in  the  TL 
curves.  In  this  section,  we  analyzed  the  curve  of  a  standard  TL  material  (LiF:Mg,Ti)  using  a 
superposition  of  first-order  TL  peaks  (Randall  and  Wilkins,  1945)  and  investigated  whether  the 
first-order  model  is  capable  of  reproducing  the  effect  of  temperature  to  which  the  material  was 
exposed  on  the  TL  curves. 

We  are  particularly  interested  in  this  study  in  using  a  model  obtained  from  TL  data  recorded  at 
low  heating  rate  and  apply  the  model  to  TL  data  recorded  at  a  higher  heating  rate.  The  goal  is  to 
extrapolate  this  approach  in  the  future  for  the  extremely  fast  heating  rates  experienced  during 
explosive  events. 


4.2. 1 . 1  The  first-order  TL  model  for  multiple  peaks 

The  first-order  TL  model  assumes  that  trapping  centers  do  not  interact,  i.e.,  the  trapped  electron 
escaping  the  trap  recombines  immediately  with  the  hole  resulting  in  luminescence.  For  a  single 
type  of  defect  characterized  by  an  activation  energy  E  and  frequency  factor  s,  the  rate  of  change 
in  the  concentration  of  trapped  charges  n  is  given  by  (McKeever,  1985): 


dn  —F/kT 

—  =  -me  '  . 
dt 


(1) 


This  equation  can  be  integrated  to  obtain  the  number  of  trapped  charges  at  any  instant  t  when  the 
material  is  subjected  to  a  temperature  profile: 

'  a 

-sje  kT(,;>ai  .  (2) 


n(t)  =  n0  exp 


In  this  equation,  «o  is  the  initial  concentration  of  charges  in  the  trapping  centers. 

During  a  TL  readout,  the  TL  intensity  (s"1)  is  proportional  to  the  rate  of  recombination.  In  quasi¬ 
equilibrium  (charges  escaping  the  trap  recombine  immediately),  the  TL  intensity  can  be  written 


as: 


■^77.(0  — 


dm 

dn 

dt 

dt 

=  me 


-F./kT 


(3) 


For  a  linear  heating  rate  such  that  T(t )  =  T0+  fit ,  Eq.  (2)  can  be  written  as 


n(T)  =  n0  exp 


(4) 


The  TL  intensity  in  K  1  then  becomes: 


P 


exp 


— —  [  exp  — — 
fi\  l  kO 


d6 


(5) 
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This  is  the  well-known  first-order  TL  model,  first  proposed  by  Randall  and  Wilkins  (1945). 

Some  important  properties  of  the  first-order  TL  processes  are  the  following: 

(a)  In  presence  of  multiple  first-order  processes,  each  process  can  be  treated  independently 
and  the  total  TL  intensity  is  a  superposition  of  the  TL  intensity  due  to  each  process; 

(b)  The  peak  temperature  does  not  depend  on  the  initial  trap  occupancy,  i.e.,  the  absorbed 
dose  in  the  crystal. 

Condition  (a)  allows  us  to  write  the  following  equation  for  a  material  containing  multiple  TL 
peaks,  the  z'-th  peak  characterized  by  E,  and  s,: 


i  P 


kT 


exp 


-JW-S) 


\dd 


(6) 


4.2. 1 .2  A  formalism  to  describe  the  effect  of  a  heating  profile 

In  this  study  we  assume  that,  previously  to  the  TL  readout,  the  material  is  subjected  to  a  heating 
profile  in  which  the  temperature  is  elevated  from  an  initial  temperature  To  until  a  maximum 
temperature  Tm  at  a  constant  rate  a  (K  s"1),  and  then  decreases  exponentially  with  a  time  constant 
r  until  it  reaches  the  initial  temperature  To.  This  temperature  profile  can  be  written  as 

f  T0+at  if  t<tm 

0(0  =  1  •  (7) 

l To  +  (Tm-T0)e  T  if t>tm 


where  we  define  tm  as  the  instant  in  which  the  temperature  reaches  a  maximum. 

According  to  Eq.  (2).  the  concentration  of  trapped  charges  after  the  material  returns  to  the  initial 
temperature,  n\ ,  is 

tn,  E  oo  E 


«0  =  «0  exp 


-s\e  kT(,)dt-sje  kT(,)dt 


(8) 


The  exponential  term  can  be  factored  into  two  terms,  one  expressing  the  fraction  of  trapped 
charges  released  during  heating,  f  =  f  (E,s,Tm,a) ,  and  one  expressing  the  trapped  charges 

released  during  cooling,  g  =  g(E,s,Tm,r) .  Using  the  temperature  profile  given  in  Eq.  (6),  this 
can  be  written  as 

«o  =  *of(E,  s,  Tm  ,a)g(E,s,Tm,  r ) , 


where 


f(E,s,Tm,a)  =  ex p 


Tm  B 

s  r  —e 


a  t 


dd 
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g(E,s,Tm,r)  =  exp 
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Therefore,  the  effect  of  the  temperature  profile  is  to  reduce  the  original  concentration  of  trapped 
charges  by  a  factor /during  the  heating  stage  times  a  factor  g  during  the  cooling  stage. 


For  a  material  containing  a  superposition  of  first-order  TL  peaks,  the  effect  of  the  temperature 
profile  to  which  the  material  was  exposed  prior  to  the  TL  readout  is  to  reduce  each  peak  by  a 
factor  figh  where  the  sub-index  /  refers  to  a  particular  TL  peak,  i.e.,  to  a  particular  set  of  E,  and  s, 
values.  The  total  TL  intensity  is  then  given  by 


HT)  ~  X  f&t  ~n~e 


exp 


-if. 

Pi 


exp|  -fe 


dd 


(12) 


4.2. 1 .3  Experimental  tests 


Experimental  TL  data  was  obtained  using  LiF:Mg,Ti  dosimeters  of  dimensions  3  mm  x  3  mm  x 
0.9  mm.  A  4-mm  thick  Kopp  5850  filter  and  a  3-mm  thick  HA3  filter  were  used  in  front  of  the 
PMT  to  reduce  the  blackbody  radiation  and  transmit  the  main  luminescence  band  from 
LiF:Mg,Ti  at  410-415  nm  (McKeever  et  al.,  1995).  No  annealing  procedure  was  used  in  between 
experiments. 


Deconvolutions  of  the  TL  curves  into  first-order  TL  peaks  were  carried  out  using  PeakFit™ 
version  4  (Jandel  Scientific  Software)  with  a  user-defined  Randall- Wilkins  model  as  in  Eq.  (5), 
the  integral  being  approximated  by  (Chen  and  McKeever,  1997): 


T  E 


—  kT 2 

k0  dd  =  ——e 


JL( 

kT 


2rr2  \ 


,  2 kT  6kzT 

1 - +  — — 


03) 


Numerical  modeling  was  carried  out  using  Mathematica  5.2  (Wolfram  Research  Inc.). 


4.2. 1 .4  Analysis  of  the  LiF  :Mg,Ti  TL  curve 

Figure  21  shows  the  TL  curve  of  LiF:Mg,Ti  recorded  at  a  heating  rate  of  0.2  K/s.  The  curve  was 
fitted  with  first-order  TL  peaks,  resulting  in  the  parameters  shown  in  Table  3. 
As  one  can  observe  in  Figure  21,  the  TL  curve  can  be  well-represented  by  a  superposition  of 
first-order  TL  peaks  in  a  first  approximation.  The  parameters  obtained  shown  in  Table  3  can  now 
be  used  to  predict  the  effect  of  an  arbitrary  temperature  profile  on  the  TL  curves  of  this  material 
using  the  formalism  previously  presented  in  this  section.  Here  we  restricted  our  investigations  to 
the  temperature  profiles  described  by  Eq.  (7). 
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Temperature  (K) 

Figure  21.  TL  curve  of  LiF:Mg,Ti  fitted  with  a  superposition  of  first-order  TL  peaks  (heating  rate  of  0.2  K/s; 
see  fitted  parameters  in  Table  3). 


Table  3.  Parameters  obtained  by  fitting  the  TL  curve  of  LiF:Mg,Ti  with  first-order  TL  peaks. 


TL  peak 

n„  (arb.  units) 

E(eV) 

l  Wl  UV1  JL  JU 

s  (s'1) 

Peak  1 

1.19x1 06 

0.740 

1.19x10'° 

Peak  2 

4.60*  106 

1.042 

2.91  xlO12 

Peak  3 

5.1  lxlO6 

1.182 

6.68xl012 

Peak  4 

5.15xl06 

1.698 

5.64xl017 

Peak  5 

1.44xl07 

1.962 

2.61xl019 

4.2.1 .5  Comparison  of  experimental  and  numerical  simulations  of  TL  curves 

Figure  22  shows  a  comparison  of  experimental  and  simulated  TL  curves  obtained  immediately 
after  irradiation  or  after  irradiation  and  pre-heating  to  different  temperatures  (373  K,  423  K,  and 
473  K)  at  a  heating  rate  of  5  K/s.  Note  that  the  heating  rate  used  in  the  pre-heating  and 
simulation  is  ~20  times  faster  than  the  heating  rate  used  to  record  the  data  from  which  the  model 
was  based  on  (0.2  K/s).  Cooling  is  assumed  to  be  immediate  in  the  case  of  both  experiment  and 
numerical  simulations. 

First  of  all,  the  results  in  Figure  22  shows  the  effect  of  the  different  effect  of  the  pre-heating  on 
the  various  TL  peaks  of  LiF:Mg,Ti.  The  important  point,  however,  is  to  note  that  the  qualitative 
agreement  between  the  model  prediction  (column  on  the  right)  with  the  experimental  data 
(columns  on  the  left).  The  implications  is  that,  in  spite  of  the  oversimplified  assumption  of  first- 
order  behavior,  the  TL  model  based  on  data  recorded  at  slow  heating  rate  was  able  to  predict  the 
effect  of  a  fast  heating  profile.  It  remains  to  be  demonstrated  whether  this  model  also  applies  to 
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the  extremely  fast  heating  profiles  experienced  during  explosive  events  and  reproducible  using 
the  microheaters  developed  in  this  project. 


Temperature  (K) 


Figure  22.  Experimental  and  simulated  TL  curves  of  LiF:Mg,Ti  (heating  rate  of  0.2  K/s)  immediately  after 
irradiation  or  subjected  to  a  pre-heating  to  different  temperatures  (373  K,  423  K,  and  473  K)  at  a  fast  heating 
rate  5  K/s.  The  simulated  curves  are  prediction  of  the  effect  of  heating  rate  on  the  TL  curves  based  on  the 
first-order  model  parameters  presented  in  Table  3. 
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4.2. 1 .6  Effect  of  different  temperature  profiles 

One  question  that  can  be  answered  based  on  this  model  is  if  a  superposition  of  TL  peaks  can 
distinguish  between  two  different  heating  profiles  that  would  affect  one  of  the  TL  peaks  in  the 
same  way.  Consider  for  example  a  heating  profile  with  maximum  temperature  of  443  K  and 
heating  rate  of  2  K/s.  Based  on  the  data  on  Table  3  and  the  formalism  developed  earlier,  such 
heating  profile  reduces  TL  peak  3  to  20%  of  its  original  value.  A  heating  profile  with  maximum 
temperature  of  406  K  and  heating  rate  of  0. 1  K/s  would  also  reduce  TL  peak  3  to  20%  of  its 
original  value.  This  means  that  in  a  material  with  a  single  TL  peak  similar  to  peak  3,  these  two 
heating  profiles  could  not  be  distinguished  based  on  the  TL  curve. 

However,  with  a  superposition  of  TL  peaks  such  as  the  one  found  in  LiF:Mg,Ti,  the  numerical 
simulations  shown  in  Figure  23a  indicates  that,  although  the  effect  of  the  two  temperature 
profiles  on  peak  3  is  identical,  the  effect  on  peak  4  is  not.  As  a  result,  a  difference  in  the  TL 
curves  can  be  observed  around  443  K.  Moreover,  this  prediction  was  confirmed  by  experimental 
data  presented  in  Figure  23b,  in  which  TL  curves  of  LiF:Mg,Ti  subjected  to  these  two 
temperature  profiles  are  compared  with  the  TL  peak  when  no  pre-heating  is  applied.  Again,  a 
difference  in  the  TL  intensity  can  be  observed  around  443  K. 


Figure  23.  (a)  Numerical  simulation  and  (b)  experimental  data  on  the  effect  of  two  heating  profiles  (maximum 
temperature  of  443  K  at  2  K/s,  and  maximum  temperature  of  406  K  at  0.1  K/s). 


4.2. 1 .7  Influence  of  cooling  stage 

The  discussion  in  the  previous  section  does  not  include  the  effect  of  the  cooling  state.  As  one 
may  expect,  the  cooling  state  will  be  important  if  cooling  is  slower  than  the  heating  stage.  Here 
we  establish  a  more  precisely  condition  to  whether  the  cooling  state  will  or  will  not  affect  the 
temperature  determination. 


Comparing  Eq.  (10)  and  (11)  one  can  see  that  the  condition  for  the  heating  stage  to  be 
predominant  is: 


(14) 
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where  T  is  in  the  range  of  temperatures  in  which  we  are  interested.  For  example,  for  a  slow 
heating  rate  of  a=  1  K/s  and  a  slow  cooling  rate  with  time  constant  r=  10  s,  the  depletion  is 
going  to  be  dominated  by  the  heating  stage.  For  a=  1  K/s  and  T ~  450  K,  this  conditions  reduces 
to  r  «  450  s.  For  a  fast  heating  rate  of  a=  103  K/s  and  T  ~  450  K,  this  conditions  reduces  to  r 
«  4.5  s. 

The  model  developed  in  the  previous  section  can  also  be  used  to  investigate  the  effect  of  cooling 
rate  in  the  TL  peaks.  Figure  24  shows  simulated  TL  curves  of  LiF:Mg,Ti  heated  to  443  K  at  a 
heating  rate  of  2  K/s,  but  in  which  the  cooling  constant  r  assumed  different  values.  According  to 
the  condition  expressed  by  Eq.  (13),  the  heating  stage  ceases  to  be  dominant  only  when  r~ 
220  s.  This  is  confirmed  by  the  simulations  presented  in  Figure  24,  which  shows  that  the  cooling 
constant  have  no  effect  on  the  TL  curves  until  it  reaches  a  value  -300  s. 


Figure  24.  Numerical  stimulation  on  the  effect  of  cooling  constant  ron  the  TL  curves  of  LiF:Mg,Ti  heated  to 
443  K  at  a  heating  rate  of  2  K/s.  The  inset  shows  the  various  simulated  temperature  profiles  to  which  the 
material  was  subjected. 


4.2.2  Implementation  of  automated  initial  rise  method 

The  model  above  demonstrates  that  it  is  possible  to  predict  the  effect  of  a  fast  heating  profile 
using  TL  parameters  (activation  energy  and  frequency  factor)  obtained  using  a  slow  heating 
profile  in  laboratory.  To  obtain  these  parameters,  the  TL  curve  must  be  fitted  with  a  set  of  TL 
peaks.  One  of  the  useful  techniques  to  obtain  the  activation  energies  of  a  TL  curve  is  to  use  a 
step-annealing  technique  and  apply  the  initial  rise  method  (McKeever,  1985).  In  this  section  we 
tested  an  automated  analysis  to  apply  the  initial  rise  method. 

The  automated  initial  rise  method  analysis  was  implemented  in  Mathematica  5.2  and  applied  to  a 
few  materials.  Figure  25  shows  representative  examples  of  TL  curves  obtained  using  partial 
cleaning  (step-annealing)  of  the  TL  curves,  which  consists  in  heating  to  the  specified  temperature 
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to  clean  the  low  temperature  TL  peaks,  followed  by  cooling  to  room  temperature  and 
measurement  of  the  whole  TL  curve.  This  procedure  subsequently  empties  the  traps  related  to 
low  temperature  peaks,  “exposing”  the  high  temperature  peaks.  The  application  of  the  initial  rise 
method  to  these  curves  allows  one  to  find  the  activation  energy  for  the  trapping  centers 
associated  with  peaks  at  increasing  temperatures,  therefore  providing  an  “activation  energy 
spectrum”  for  the  TL  in  the  material.  One  can  see  in  Figure  25  that  the  automated  method 
produced  activation  energies  that  are  very  similar  to  the  manual  analysis  of  the  curves,  although 
some  discrepancies  were  observed  for  the  YO  sample  (Figure  25b). 

In  order  to  test  the  reliability  of  this  analysis  procedure,  we  also  simulated  a  material  containing 
four  first-order  TL  peaks  (Figure  26).  The  simulated  data  was  then  treated  as  the  experimental 
data  and  analyzed  using  the  initial  rise  method.  The  results  showed  in  Figure  27  demonstrate  that 
under  these  ideal  conditions  the  activation  energies  could  be  successfully  recovered.  However,  in 
regions  where  the  TL  curves  have  low  intensity,  anomalous  values  of  activation  energies  were 
obtained. 

The  conclusion  is  that  an  ideal  situation  is  to  materials  with  well-defined  first-order  TL  peaks. 
This  would  make  it  easier  to  obtain  the  TL  parameters  and  done  for  LiF:Mg,Ti  in  Section  4.4.2. 
If  not  possible,  the  approach  described  here  could  be  used  to  extract  the  parameters  for  a  large 
superposition  of  TL  peaks. 


Figure  25.  Application  of  initial  rise  method  to  Gd2SiOs:Ce  and  Y203  (undoped). 
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Figure  26.  Simulated  TL  curve  containing  four  peaks  with  activation  energy  of  0.8  eV,  0.9  eV,  1.2  eV  and 
1.2  eV. 


Simulation:  4  peaks  E  ='{0.8eV,  0.9eV,  1.1  eV,  1.2eV} 


Figure  27.  Activation  energy  recovered  using  the  automated  initial  rise  method  and  manual  method. 


4.3  Experimental  Demonstration  of  Microparticle  Thermal  Sensors 


4.3.1  Background 

The  investigation  and  demonstration  of  temperature  sensing  using  thermoluminescent 
microparticles  was  primarily  performed  at  the  University  of  Minnesota  by  Ph.D.  graduate  student 
Merlin  Mah  under  the  direction  of  Professor  Talghader  with  support  from  Professor  Yukihara  of 
OSU. 
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Thermoluminescent  materials  have  extensive  applications  in  radiation  dosimetry 
(Escobar-Alarcon,  2003)  and  sedimentary  dating  (Wintle,  1979).  In  dosimetry,  ionizing  radiation 
excites  electron-hole  pairs  in  high  bandgap  materials,  filling  previously-empty  electron  traps; 
thermoluminescence  (TL)  is  later  used  to  assess  their  occupancy.  In  dating,  background  radiation 
from  sediment  fills  the  traps  over  time  (the  sun  having  emptied  all  traps  prior  to  burial). 

The  lifetime  of  trapped  charges  is  determined  by  the  Arrhenius  equation,  which  states  the 
probability  of  detrapping  per  unit  time  as  equal  to  s  exp  (- E/kT ),  where  s  is  a  frequency  factor 
measured  in  s'1,  E  is  the  activation  energy  for  the  process,  T  is  the  temperature  and  k  is  the 
Boltzmann  constant  (McKeever,  1985).  Since  the  Arrhenius  relationship  includes  a  strong 
variation  with  temperature,  the  luminescence  of  a  material  should  also  serve  as  a  probe  of  the 
thermal  history  to  which  that  material  has  been  exposed.  This  concept  has  only  previously  been 
applied  to  natural  TL  materials  in  meteorites  to  determine  their  probable  orbital  distance  prior  to 
entry  into  Earth’s  atmosphere  (McKeever,  1980)(Benoit,  1991)(Sears,  1974). 

It  is  demonstrated  in  this  paper  that  thermoluminescent  materials  can  be  used  as  active 
temperature  sensors.  Microparticles  of  TL  materials  are  uniquely  suited  for  the  detailed  profiling 
of  extremely  violent  temperature  events,  such  as  the  explosive  destruction  of  biological  particles, 
that  would  destroy  existing  direct-contact  sensors.  TL  microparticles  have  no  fabricated 
components,  are  often  composed  of  highly  robust  materials,  are  inexpensive  and  relatively 
simple  to  fabricate  in  large  quantities,  can  be  micromachined  to  match  the  size  and  shape  of  a 
bioagent  or  to  be  injected  into  an  otherwise-sealed  system,  and  can  be  spread  en  masse  in  order 
to  probe  temperature  differences  in  different  spatial  regions. 


4.3.2  Microparticles  and  Microheaters 

Mg2SiC>4:Tb,Co  (magnesium  silicate)  microparticles  were  chosen  for  having  two  or  more 
well-separated  thermoluminescence  peaks  occurring  at  conveniently  reachable  temperatures. 
Microparticles  with  diameters  on  the  order  of  10-20  pm  were  prepared  (Mittani,  2008)  by  Dr. 
M.  Prokic,  Institute  of  Nuclear  Sciences.  They  are  sprinkled  onto  microheaters,  as  shown  in 
Figure  28. 


HDTRA1  -07-1  -001 6 

University  of  Minnesota  (Twin  Cities) 


(a)  A  microheater  topped  with  several  Mg2Si04:Tb,Co  microparticles 


(b)  A  close-up  of  the  Mg2Si04:Tb,Co  microparticles 

Figure  28  -  Several  Mg2Si04:Tb,Co  microparticles  on  a  microheater  platform,  with  many  more  scattered 
nearby. 
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Microheaters  are  used  because  they  are  able  to  rapidly  and  repeatedly  bring  the 
microparticles  from  room  temperature  to  over  500  °C  and  back  with  time  constants  on  the  order 
of  tens  of  milliseconds — 32  ms  in  our  current  devices,  as  measured  by  bolometric  methods 
(Kruse,  1962)  —  as  well  as  provide  a  standard  linear  TL  ramp.  They  are  composed  of  a  200  nm 
thick  etch-released,  corner-supported  platform  of  LPCVD  silicon  nitride  topped  with  a  1 60  nm 
Ti/Pt  serpentine  resistor.  The  structure  designs  vary  considerably  in  size  but  the  data  reported 
here  use  heaters  with  platforms  of  side  length  300  pm. 

The  platinum  exhibits  a  highly  linear  resistance-to-temperature  relationship,  which  is 
extrapolated  by  heating  the  entire  device  and  substrate  on  a  hotplate  while  measuring  the 
resistance;  thereafter,  the  operating  temperature  of  the  microheater  can  be  monitored  and 
controlled  via  the  resistance  of  the  platinum,  just  as  in  a  common  thermistor.  Figure  29  shows  a 
set  of  such  measurements  for  a  microheater  on  the  die  used  here.  The  microheater  is  assumed  to 
be  at  the  same  temperature  as  the  hotplate  surface,  measured  by  a  R.T.D.  Co.  2823  resistive 
temperature  device  located  adjacent  to  the  microheater  die  on  the  hotplate.  Only  the  temperature 
range  from  35  °C  to  60  °C  is  considered:  data  below  this  range  are  omitted  due  to  insufficiently 
significant  heating,  and  data  above  it  are  excluded  because  of  convection  effects.  The  maximum 
deviation  among  the  independent  linear  fits  of  the  data  sets  when  evaluated  at  350  °C  is  15.57 
°C,  or  4.4%.  The  microheater  resistance  is  measured  using  a  Keithley  2410  sourcemeter,  which 
will  also  serve  to  drive  and  monitor  the  microheater  in  all  pulsed  and  linear  ramp  heating 
reported  herein. 


Hotplate  temperature  (°C) 

Figure  29  -  Three  sets  of  data  and  the  resulting  average  linear  fit  for  the  resistance  of  a  hotplate-heated 
microheater.  Each  point  represents  the  average  of  5  measurements  taken  in  rapid  succession;  no  other 
noise  reduction  has  been  used.  To  reduce  the  effects  of  inconstant  hotplate  heating,  the  measurements  are 
taken  as  the  powered-down  hotplate  cools  to  ambient  temperature. 
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Prior  to  testing,  the  traps  of  the  microparticles  are  filled  by  irradiating  with  UV  light  in 
situ  on  the  microheater  by  exposure  to  a  broad-spectrum  1 80-400  nm  30  W  deuterium  lamp.  A 
200  ms  pulse  of  variable  peak  temperature,  measured  using  the  platinum  thermistor,  is 
administered  to  a  magnesium  silicate  particle  with  filled  traps.  The  thermoluminescence  glow 
curve  is  then  gathered  using  a  standard  (McKeever,  1985)(Furetta,  2003)  linear  temperature 
ramp  of  2  °C/sec,  reaching  to  350  °C.  The  emitted  photons  are  counted  using  a  Sens-Tech 
P25PC  photomultiplier  tube  interfaced  to  a  BK  1 823A  frequency  counter.  The  entire  process  is 
automated  to  run  in  a  room  darkened  such  that  the  PMT  background  count  is  three  times  its 
specified  maximum  dark  count.  This  also  eliminates  any  effects  due  to  photobleaching  that  are 
experienced  by  many,  but  not  all,  1L  materials  (McKeever,  1985).  For  display,  each  raw  data 
point  is  averaged  with  its  20  closest  neighbors  (corresponding  to  approximately  1 0  °C)  and  each 
curve  multiplicatively  scaled  so  that  their  heights  are  equal  at  266  °C,  the  apex  of  the  second 
peak.  The  resulting  data  evinces  a  strong  background  intensity,  likely  from  thermal  blackbody 
emission  and  nearby  TL  signatures,  which  becomes  dominant  at  ramp  temperatures  above  310 
°C  after  the  Mg2Si(>4:Tb,Co  TL  peaks  have  died  away.  This  unwanted  signal  is  modeled  as  an 
exponential  function  curvefitted  to  the  non-pulsed  Mg2Si04:Tb,Co  above  310  °C,  and  subtracted 
from  the  intensity  data  of  each  subsequent  pulsed  run.  The  actual  TL  signatures  of  interest  can 
then  be  examined  directly.  Figure  30  compares  the  curve-matched  background  function  to  the 
intensity  data  before  and  after  this  background  is  subtracted. 


Figure  30  -  The  intensity  observed  during  the  collection  of  the  TL  glow  curves  before  and  after 
subtraction  of  curve-fitted  exponential  background  function. 
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The  resulting  thermoluminescence  intensity  of  Mg2Si04:Tb,Co  versus  linear  ramp  temperature 
for  a  series  of  different  thermal  histories,  in  this  case  peak  pulse  temperatures,  is  shown  in  Figure 


Heater  temperature  CC) 


Figure  31  -  The  thermoluminescence  glow  curves  of  Mg2Si04:Tb,Co  microparticles  after  a  °190  ms 
explosive  heating  pulse.  The  legend  at  right  indicates  the  peak  temperature  of  the  pulse  corresponding  to 
each  curve.  The  intensities  shown  are  normalized  at  266  °C,  the  top  of  the  second  observable  peak.  The 
ratios  of  the  intensities  of  the  peaks  strongly  correlate  with  the  maximum  temperature  to  which  the 
particles  had  been  exposed. 


4.3.4  Discussion  and  Theory 


The  exact  mechanics  and  trap  parameters  of  Mg2Si04:Tb,Co  are  unknown,  and  indeed 
have  been  observed  to  vary  greatly  according  to  preparation.  We  therefore  modeled  the  material 
as  two  superimposed  traps  exhibiting  first-order  kinetics.  Under  this  model,  the  population  n  of 
an  electron  trap  of  depth  E,  frequency  factor  s,  and  initial  population  no  after  experiencing 
temperature  profile  T(t)  is  given  (Randall,  1945)  by: 

n  (t)  =  n0  exp  f—  f  se~E^kT^  dt] 

\  Jo  J  (1) 

If  this  population  value  is  taken  to  be  n0  ‘  after  T(t)  has  run  its  course,  then  the  TL  intensity 
observed  from  the  remaining  population  during  a  linear  temperature  ramp  T’(t  ’)  is  expressed  as 


^Cno  exp 


se 


—E/kT(t) 


se-E/kT'{t') 


(2) 
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where  C  is  a  constant.  Lacking  knowledge  of  the  trap  parameters  of  this  material,  we  chose  E,  s, 
and  no  values  that  best  fit  our  non-pulsed  experimental  curve;  numerical  simulations  based  on  the 
preceding  equations  were  then  used  to  produce  pulsed  glow  curves  for  these  theoretical  particles, 
figure  32  compares  the  relation  between  the  intensity  peak  height  and  temperature  pulse  maxima 
of  the  resulting  simulations  to  that  of  the  experimental  data  seen  in  Figure  31. 


Figure  32  -  The  ratio  of  the  height  of  the  first  TL  peak  to  the  height  of  the  second  as  a  function  of  pulse 
temperature  for  simulated  and  experimental  data. 


A  one-to-one  decrease  occurs  in  the  intensity  ratio  of  the  first  peak  to  second  peak 
heights  as  the  pulse  temperature  climbs.  The  slope  of  the  curve  is  reasonably  stable,  indicating 
that  the  sensitivity  of  the  material  will  be  fairly  steady  over  the  usable  pulse  temperature  range. 
In  the  range  we  observed,  model  and  data  differed  by  an  average  of  4.4%,  with  no  errors  greater 
than  9.1%.  It  is  therefore  possible  to  deduce  the  pulse  temperature  from  the  TL  glow  curve  by 
comparing  the  ratios  of  intensity  peak  heights,  provided  the  length  of  the  pulse  is  known.  Even 
better  accuracy  would  be  obtained  by  simply  curve  fitting  the  I  L  peak  ratio  in  much  the  same 
manner  that  the  resistance  of  commercial  thermistors  is  calibrated.  With  this  particular  material 
and  temperature  pulse  duration  the  first  TL  peak  becomes  indistinguishable  after  the  pulse 
temperature  passes  290  °C,  but  other  TL  materials  with  deeper  traps  could  be  substituted  or  used 
in  parallel  to  take  over  from  Mg2Si04:Tb,Co  at  this  point. 

lo  conclude  this  section,  we  have  demonstrated  the  recovery  to  good  accuracy  of  the 
peak  temperatures,  ranging  from  232  °C  to  313  °C,  experienced  in  a  microheater-powered  rapid 
heating  event  of  known  duration  less  than  200  ms.  The  observed  change  in  the 
thermoluminescence  peak  intensity  ratio  of  two  Mg2Si04:Tb,Co  traps  has  been  predicted  with  an 
average  error  of  4.4%  using  first-order  kinetics  modeling. 
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4.4  Characterization  of  existing  dosimetric  materials  for  thermometry 
4.4.1  Overview  of  the  material  research  and  literature  research 

In  spite  of  the  relatively  large  availability  of  TL  materials  for  dosimetry  application,  most  of  the 
materials  available  are  not  suitable  for  temperature  sensing.  This  stems  from  the  fact  that  the 
requirements  for  dosimetry  are  different  than  the  requirements  for  temperature  sensing.  Detectors 
that  would  be  ideal  for  temperature  sensing  were  usually  rejected  as  dosimetric  materials  in  the 
past. 

The  characteristics  of  materials  of  interest  for  temperature  sensing  are  the  following: 

(i)  High  thermoluminescence  output  (i.e.,  high  sensitivity); 

(ii)  Presence  of  several  thermoluminescent  peaks  covering  a  wide  range  of  temperatures; 

(iii)  First-order  kinetic  behavior  or  thermoluminescence; 

(iv)  Thermal  stability  of  the  thermoluminescent  peaks; 

(v)  Optical  stability  of  the  thermoluminescent  peaks  (i.e.,  TL  is  insensitive  to  light); 

(vi)  Chemical  and  mechanical  stabilility. 

There  is  currently  no  material  satisfying  all  the  criteria  above  and,  therefore,  research  into  new 
luminescent  materials  that  may  have  been  previously  ignored  for  dosimetry  applications  is  of 
paramount  importance. 

The  material  investigations  in  this  project  focused  on  silicates  and  oxides  due  to  their  chemical 
stability.  They  included  new  nanophosphors  produced  by  Luiz  B.  Jacobsohn  and  Michael  Blair 
at  Los  Alamos  National  Laboratory  (LANL),  and  natural  materials  known  to  exhibit  potentially 
good  thermoluminescence  properties  (natural  topaz).  The  main  findings  are  described  below. 


4.4.2  Commercial  dosimetric  materials 
4.4.2. 1  A1203:C 

A1203:C  is  one  of  the  most  sensitive  TL  and  OSL  material  available  with  relatively  well-known 
TL  properties  (McKeever  et  al.,  1995).  Furthermore,  its  hardness  and  resistance  to  chemical 
attack  makes  this  material  ideal  for  harsh  environments.  The  problem  is  that  the  TL  of  A1203:C  is 
dominated  by  a  peak  at  180°C,  which  is  the  peak  used  in  dosimetry.  Therefore,  to  extend  the 
application  of  the  material  to  higher  temperatures,  it  is  necessary  to  investigate  its  high- 
temperature  peaks. 

Figure  33  shows  the  TL  curves  of  A1203:C  after  irradiation  with  a  high  dose  of  300  Gy  and  pre¬ 
heating  to  280°C  to  eliminate  the  main  dosimetric  peak.  The  residual  TL  curve  shows  various 
peaks  ranging  from  300°C  to  600°C  which  can  be  useful  for  thermal  sensing. 

The  problem  of  using  TL  of  A1203:C  is  that  the  thermal  quenching,  i.e.,  reduction  in 
luminescence  efficiency  with  increasing  temperatures,  decreases  the  intensity  of  the  TL  peaks 
significantly.  The  solution  to  this  problem  is  to  try  using  optical  stimulation  to  access  these  deep 
traps,  taking  advantage  of  the  fact  that  OSL  is  measured  at  room  temperature  and,  therefore,  it  is 
not  affected  by  thermal  quenching. 
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The  possibility  of  probing  the  deep  traps  associated  with  the  high-temperature  TL  peaks  in 
A1203:C  was  investigated  using  blue  light.  Figure  34  shows  the  effect  of  illumination  on  the  TL 
peaks.  It  is  clear  from  this  data  that  blue  stimulation  can  effectively  deplete  the  charge  carriers  in 
these  deep  traps  and,  therefore,  can  be  used  to  probe  the  signal  and  obtain  a  calibration  curve 
similar  to  the  one  shown  in  Figure  43b,  but  for  the  OSL  signal  instead  of  TL  signal. 


AI203,  sample#2:  Signal  from  TL700°C  after  TL280°C ,  irradiated  300Gy 


Figure  33.  AI203:C  sample  irradiated  with  300  Gy,  and  pre-heated  to  280°C  to  remove  the  main  dosimetric 
trap  at  180°C,  revealing  the  deep  traps  present  in  the  material. 
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TL  700°C:  Contribution  of  Blue  OSL  (1800s)  to  the  measurement 


Temperature  (°C) 


Figure  34.  Effect  of  blue  stimulation  on  deep  traps  of  an  A^OjtC  single  crystal.  The  sample  was  irradiated 
with  300  Gy  and  pre-heated  to  280°C  to  remove  the  main  dosimetric  trap.  The  graph  shows  the  remaining  TL 
curve  before  and  after  blue  stimulation  for  1800  s  (LEDs  centered  at  470  nm,  ~50  mW/cm2). 


4A.2.2  LiF:Mg,Ti 

LiF:Mg,Ti  is  a  TL  material  widely  used  in  TL  dosimetry  (McKeever  et  al.,  1995).  However, 
there  are  several  properties  that  make  this  material  of  interest  for  temperature  sensing,  including 
the  presence  of  several  TL  peaks,  the  thermal  stability,  the  insensitivity  to  light  exposure,  and  the 
ready  commercial  availability.  The  disadvantage  is  the  complex  defects  associated  with  the  TL 
process,  which  may  affected  our  ability  to  extract  information  from  the  TL  curves.  Here  we 
present  some  of  the  tests  performed  with  this  material  for  application  in  temperature  sensing. 
More  detailed  investigations  were  carried  out  using  this  material  in  Section  4.4.2. 

Figure  35  shows  that  the  TL  curves  of  LiF:Mg,Ti  are  very  reproducible,  which  is  a  requirement 
for  use  of  the  material  as  a  temperature  sensor.  The  material  also  shows  high  sensitivity.  The  TL 
intensity  of  ~10  grains  reaches  an  intensity  close  to  saturation  of  the  PMT  (Figure  36).  The  TL 
peaks  above  100°C  are  stable  at  room  temperature  and  could  be  used  as  temperature  sensors 
(Figure  37). 
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Figure  35.  Reproducibility  test  of  Lit  :Mg,Ti.  The  data  shows  the  TL  curves  for  five  cycles  of  irradiation  and 
TL  readout.  The  results  demonstrate  the  reproducibility  of  the  TL  curves  from  this  material. 


Figure  36.  I  L  curves  of  grains  of  LiF:Mg,Ti  irradiated  with  a  large  dose.  This  curves  provide  an  idea  of  the 
TL  intensity  that  could  be  obtained  for  single  grains  of  the  material  to  be  used  in  temperature  sensing.  The 
signal  intensity  (~106  cps)  is  close  to  the  saturation  of  the  PMT. 
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Figure  37.  TL  curves  of  LiF:Mg,Ti  immediately  following  irradiation  and  after  1  or  2  h  after  irradiation.  The 
data  demonstrates  that  the  TL  peaks  at  temperatures  higher  than  100°C  are  relatively  stable  and  could  be 
used  to  temperature  sensing.  The  peaks  below  100°C  decay  at  room  temperature  and  are  not  suitable  for 
temperature  sensing. 


4.4.3  Materials  under  development 
4.4.3. 1  Silicates  and  aluminates 

Most  of  the  FL  materials  available  today  were  developed  for  radiation  dosimetry  applications 
and  optimized  to  have  a  single  or  easily  distinguished  TL  peak.  On  the  other  hand,  for  the 
applications  intended  in  this  project  new  materials  with  multiple,  equally  intense  TL  peaks  are 
necessary.  Therefore,  it  is  necessary  to  research  new  materials  that  could  be  used  as  “model” 
materials  to  test  our  hypothesis  and  eventually  for  the  real  application  in  temperature  sensing. 

Figure  38  shows  the  TL  curves  of  several  silicates  and  aluminates  that  we  obtained  through  other 
collaborations  and  investigated  for  potential  application  as  temperature  sensors.  The  materials 
have  very  different  distribution  of  trapping  centers,  but  the  Mg2Si04  samples  (Figure  38b)  shows 
an  interesting  combination  of  high  intensity  1 L  peaks.  Sample  16  is  particularly  interesting  due 
to  the  presence  of  at  least  three  TL  peaks  in  the  range  from  150°C  to  400°C. 

The  reproducibility  of  the  TL  curves  from  Mg2Si04:Tb,Co  (sample  16)  is  demonstrated  in  Figure 
39,  which  shows  the  TL  curves  obtained  after  five  cycles  of  readout  and  irradiation.  The 
intensity  of  the  TL  curves  do  not  deviate  by  more  than  2%. 

Figure  40  demonstrates  the  thermal  stability  of  the  TL  peaks  under  normal  conditions  of  storage 
at  room  temperature.  Except  for  the  low  temperature  portion  of  the  TL  curve  <100°C,  the  TL 
peaks  are  relatively  stable  over  up  to  24  h  of  storage. 

Figure  41  shows  the  TL  curves  of  Mg2Si04  sample  16  as  a  function  of  maximum  temperature  of 
the  pre-heating.  The  behavior  of  the  TL  peaks  is  similar  to  the  model  predictions,  showing 
decrease  of  the  TL  peaks  at  different  temperatures  (see  Figure  3).  If  the  TL  peak  intensity  is 
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plotted  against  the  maximum  pre-heating  temperature,  the  graph  shown  in  Figure  42  is  obtained. 
This  graph  allows  one  to  determine  the  maximum  temperature  reached  by  the  crystal  in 
laboratory  conditions. 

The  procedure  for  determining  the  maximum  temperature  is  illustrated  in  Figure  43.  Suppose 
that  a  material  irradiated  with  a  known  dose  is  heated  to  an  unknown  temperature  and,  after 
recovery,  the  TL  curve  is  measured  (dashed  line).  The  material  can  then  be  irradiated  again  with 
the  same  dose  to  record  the  original  TL  curve  (solid  line).  Comparison  between  these  two  curves 
(Figure  43a)  together  with  the  calibration  curve  obtained  in  laboratory  (Figure  43b)  allows  the 
maximum  temperature  to  be  determined. 

This  simple  approach  illustrates  the  usefulness  of  the  technique,  but  more  advanced 
mathematical  procedure  such  as  glow  curve  deconvolution. 

In  this  particular  case,  the  calibration  curve  is  valid  only  for  the  heating  rate  used  in  the 
experiments.  However,  as  shown  in  the  simulations,  the  variation  in  the  calibration  curve 
decreases  as  the  heating  rate  increases.  The  calibration  curve  can  then  be  obtained  for  fast 
heating  rates  and  realistic  time-temperature  profiles  using  the  microheaters  developed  in  the 
present  project. 
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Figure  38.  TL  curves  of  10  different  new  TL/OSL  materials.  The  multiple  TL  peaks  in  the  MgsSi04  samples 
makes  them  a  suitable  materials  for  testing. 


REPRODUCIBILITY  CHECK 


Figure  39.  Reproducibility  of  the  TL  curves  from  Mg2Si04:Tb,Co  (sample  16)  after  five  cycles  of  readout  and 
irradiation  with  the  same  dose. 
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Figure  40.  Short  term  stability  of  the  TL  peaks  from  sample  Mg2Si04:Tb,Co  (sample  16).  The  sample  was 
irradiated  and  TL  readout  was  performed  after  various  interval  of  time  following  irradiation. 


Figure  41.  TL  curves  of  MgSi04  samples  after  irradiation  with  300  Gy  and  pre-heating  to  increasing 
maximum  temperatures  Tmax.  Each  curve  corresponds  to  a  specific  maximum  temperature  TmtLx. 
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Figure  42.  Intensity  of  the  two  prominent  TL  peaks  in  MgSiCL  as  a  function  of  maximum  pre-heating 
temperature. 


Figure  43.  Illustration  of  one  approach  to  determine  the  maximum  temperature  experienced  by  the 
luminescence  material. 


4.4. 3.2  Oxyorthosilicates  (YSO,  GSO,  LSO) 

A  comparison  was  made  between  powder  obtained  from  materials  produced  by  single  crystal 
growth  (Czochraslki  method)  and  nanophosphors  produced  by  Solution  Combustion  Synthesis 
(SCS).  Table  4  summarizes  some  of  the  characteristics  of  the  samples,  where  the  most  important 
result  is  the  radioluminescence  (RL)  intensity.  This  measurement  demonstrates  that  the  intensity 
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of  the  radioluminescence  of  the  nanophosphors  is  comparable  to  oxyorthosilicates  produced  by 
single  crystal  growth,  indicating  that  the  Solution  Combustion  Synthesis  can  successfully 
produce  materials  of  comparable  light  yield  for  scintillation  applications. 

Figure  44  shows  the  thermoluminescence  curves  of  the  various  samples,  normalized  to  the  mass 
of  material  used  in  the  measurements.  The  curves  show  first-order  behavior,  as  indicated  by  the 
fact  that  the  curve  shape  is  independent  on  the  dose  (see  insets  in  Figure  44).  Flowever,  the 
thermoluminescence  intensity  of  the  nanophosphors  is  relatively  low  compared  to  the  bulk 
crystals.  This  indicates  a  low  concentration  of  trapping  centers. 

Another  disadvantage  is  the  lack  of  thermal  stability  of  the  TL  curves.  Figure  45  shows  the  TL 
curves  for  various  samples  for  different  times  after  irradiation.  Except  for  the  YSO  produced  by 
Czochralski  method,  all  samples  show  a  significant  fading  even  in  the  high  temperature  regions 
of  the  TL  curve  (see  Figure  45b). 

The  oxyorthosilicate  samples  also  exhibit  strong  sensitivity  to  light,  as  shown  in  Figure  46.  Light 
exposure  eliminates  most  of  the  TL  signal  induced  by  radiation.  This  characteristic  hampers  the 
use  of  these  materials  in  experiments  where  light  exposure  cannot  be  avoided. 

In  summary,  oxyorthosilicate  nanophosphors  have  shown  to  have  excellent  properties  for  use  as 
scintillating  materials  in  radiation  detection,  but  their  usefulness  for  temperature  sensing  seems 
to  be  limited  by  the  low  concentration  of  trapping  centers,  strong  fading  of  the  signal,  and  light 
sensitivity.  However,  we  should  not  exclude  the  possibility  of  finding  other  nanophosphors  that 
do  not  show  the  properties  above. 


Table  4.  Main  properties  and  characteristics  of  bulk  and  nanophosphor  samples  of  LSO,  GSO,  and  YSO  used 
in  this  study.  


Property/detail 

LSO 

(bulk) 

SCS-LSO- 

24 

GSO 

(bulk) 

SCS-GSO-8 

YSO 

(bulk) 

SCS-YSO- 

26 

Synthesis 

Czochralski 

scs 

Czochralski 

SCS 

Czochralski 

SCS 

Grain  size  range 

125-250  pm 

20-100nm 

75-125  jam 

20-100nm 

45  jam 

20-100nm 

RL  intensity 

(103  counts  s'1  mg’1)^ 

49.0  ±2.4 

64  ±4 

94  ±  13 

36  ±4 

97  ±  14 

44.3  ±  1.5 

Main  TL  peaks(b) 

340  K  (I) 

405  K  (II) 
480  K  (III) 
565  K  (IV) 

345  K 

373  K 

365  K 

398  K 

365  K 

E  (eV)(c) 

1.129 

0.782 

1.016 

1.009 

1.005 

0.987 

1.73  x  1015 

2.07  x  1010 

4.42  x  1012 

8.14  x  1012 

3.83  x  10" 

3.58  x  1012 

<a)  Radioluminescence  measured  under  X-ray  excitation  (35  kVp)  using  Schott  BG-39  filters  in  front  of  the  PMT. 

The  values  are  the  mean  and  experimental  standard  deviation  of  the  values  obtained  for  three  aliquots  with  weight 
between  1.5  and  10  mg. 

<b)  Based  on  TL  curves  measured  at  l°C/s. 

™  Activation  energy  for  peak  I  determined  using  the  Initial  Rise  Method. 

Frequency  factor  for  peak  I  calculated  based  on  the  activation  energy  and  temperature  of  the  peak,  assumin°  a 
first-order  TL  process. 
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Figure  44.  TL  curves  of  oxyorthosilicates  produced  by  the  Czochralski  method  as  compared  to 
oxyorthosilicates  produced  by  the  Solution  Combustion  Synthesis  (SCS)  method.  (See  text  for  more  details.) 
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Figure  45.  TL  curves  of  oxyorthosilicates  produced  by  the  Czochralski  method  and  nanophosphors  produced 
by  the  Solution  Combustion  Synthesis  (SCS)  method  for  various  intervals  between  irradiation  and 
measurement. 
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Figure  46.  Light  sensitivity  (blue  LED  illumination)  of  YSO  oxyorthosilicate  samples  prepared  by  (a) 
Czochralski  method  and  (b)  Solution  Combustion  Synthesis  method. 


4.4.4  Natural  materials 

An  alternative  route  to  finding  materials  with  suitable  properties  for  this  project  is  to  investigate 
natural  crystals.  The  advantage  is  that  natural  crystals  have  a  variety  of  trapping  centers  and, 
whereas  the  concentration  and  type  of  defects  may  vary,  properties  associated  with  the  matrix 
may  indicate  the  type  of  materials  that  could  be  produced  in  the  future.  Natural  topaz, 
Al2Si04(F,0H)2,  is  a  material  that  has  been  investigated  for  radiation  dosimetry  and  is  known  to 
have  TL  emission  (Yukihara  and  Okuno,  1998). 

The  samples  used  in  this  study  are  two  natural  samples  from  Brazil,  labeled  JMI01  and  JOZ07. 
The  TL  curves  of  these  samples  are  shown  in  Figure  47,  where  one  can  see  one  excellent 
characteristic  for  temperature  sensing:  the  presence  of  several  TL  peaks  in  a  temperature  range 
from  room  temperature  up  to  500°C.  The  TL  curves  exhibited  good  reproducibility,  as  shown  in 
Figure  48.  The  TL  curves  grow  with  dose  (Figure  49)  and  the  curve  shape  depends  only  slightly 
on  the  irradiation  dose  (Figure  50),  except  for  the  high  temperature  region  which  is  affected  by 
the  background  emission  due  to  blackbody  radiation.  This  indicates  that  the  curves  can  be 
described  by  a  superposition  of  first-order  TL  peaks,  as  required  in  this  project. 

The  low  temperature  TL  peaks  are  naturally  unstable  at  room  temperature,  but  Figure  51  shows 
that  the  high  temperature  peaks  are  not  affected  by  blue  light.  So,  as  opposed  to  all  materials 
previously  investigated,  topaz  seems  to  be  light  insensitive.  In  addition,  Figure  52  shows  that  TL 
can  also  be  induced  by  UV  illumination,  although  in  this  case  the  TL  curve  is  significantly 
different  from  the  beta  induced  TL  curves.  More  investigations  need  to  be  performed  to 
determine  whether  UV  can  be  used  or  not  to  fill  the  traps  in  the  intended  application. 

The  only  drawback  observed  so  far  is  the  low  light  intensity  of  the  TL  from  natural  topaz.  Table 
5  compares  the  TL  light  yield  from  topaz  as  compared  to  A1203:C,  where  one  can  see  that  topaz 
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is  at  least  three  orders  of  magnitude  less  bright  than  A1203:C.  Nevertheless,  it  should  be  pointed 
out  that  this  is  not  surprising,  since  ALOstC  is  one  of  the  most  sensitive  luminescent  dosimeters 
available  today.  As  long  as  one  can  measure  sufficient  TL  signal  from  the  grains  to  be  used  in 
the  temperature  measurements,  the  overall  sensitivity  should  not  be  a  problem. 

In  spite  of  the  low  TL  intensity,  the  natural  topaz  samples  have  shown  good  properties  to  serve 
as  a  test  materials:  reproducibility,  almost  first-order  behavior,  and  light  insensitivity.  We  intend 
to  perform  tests  for  the  intended  application  using  this  sample,  at  least  until  a  more  suitable 
material  is  found  or  produced. 


Figure  47.  TL  curves  of  two  topaz  crystal  samples,  JMI01  (single  crystal,  28mg)  and  JOZ07  (powder,  20mg). 
The  TL  measurements  were  carried  out  at  l°C/s  using  Kopp  5850  +  HA3  optical  fdters  in  front  of  the  PMT. 
The  samples  were  irradiated  with  beta  rays  from  a  "Sr/^Y  source.  The  TL  intensities  were  normalized  to  the 
mass  of  the  sample  and  dose  of  radiation. 


Figure  48.  TL  curves  of  topaz  sample  JMI01  (single  crystal,  28  mg)  repeatedly  irradiated  with  3  Gy.  The  TL 
measurements  were  carried  out  at  l°C/s  using  Kopp  5850  +  HA3  optical  filters  in  front  of  the  PMT. 
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Figure  49.  TL  curves  of  topaz  sample  JMI01  (single  crystal,  28  mg)  and  JOZ07  (powder,  20  mg)  for  various 
irradiation  doses  up  to  6  Gy.  The  TL  measurements  were  carried  out  at  l°C/s  using  Kopp  5850  +  HA3  optical 
filters  in  front  of  the  PMT. 


Figure  50.  Normalized  TL  curves  of  sample  JMI01  (single  crystal,  28  mg)  irradiated  with  various  doses  (data 
shown  in  Figure  49a). 
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Figure  51.  TL  curve  of  irradiated  topaz  samples  JMI01  (single  crystal,  28  mg)  and  JOZ07  (powder,  20  mg)  as 
a  function  of  blue  illumination  time.  The  samples  were  irradiated  with  3  Gy  of  beta  rays  from  a  90Sr/90V 
source  and  dluminated  with  blue  light  from  LEDs  (470nm,  20  mW/cm2).  TL  measurements  were  carried  out 
at  l°C/s  using  Hoya  U-340  optical  filters  in  front  of  the  PMT. 


Figure  52.  TL  curves  of  topaz  sample  JOZ07  (powder,  20  mg)  after  UV  illumination  and  beta  irradiation.  The 
UV  illumination  consisted  in  exposure  to  unfiltered  light  from  a  deuterium  lamp  (2.44  mW/cm2  at  the  sample 
position)  for  5  min,  and  the  beta  irradiation  dose  was  done  using  a  "Sr/^Y  source.  The  TL  measurements 
were  carried  out  at  l"C/s  using  Kopp  5850  +  HA3  optical  filters  in  front  of  the  PMT. 


Table  5.  Light  yield  comparison  between  topaz  samples  and  dosimetric  A12Q3:C.  The  samples 

Material/Sample  TL  yield  Relative  TL  yield 

(counts  mg'  mGy’) (%) 


AUOy.C 

Topaz/JMIOl 

Topaz/JQZ07 


1920 

2.57 

1.37 


100 

0.13 

0.07 
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4.4.4. 1  Kyanite,  sillimanite,  and  andulusite 

Kyanite,  sillimanite,  and  andalusite  are  polymorphs  of  aluminum  silicate  (AbSiOs).  They  were 
investigated  here  to  explore  the  TL  properties  of  this  type  of  material.  Unfortunately,  these 
materials  showed  only  low  temperature  TL  peak  (<1 00°C)  which  is  not  suitable  for  temperature 
sensing. 


Figure  53.  TL  curves  of  kyanite,  sillimanite,  and  andalusite,  which  are  natural  forms  of  AI2SiOs,  after 
irradiation  with  1  Gy. 


4.5  Development  of  new  materials  for  luminescence  thermometry 
4.5.1  Solution  Combustion  Synthesis 

In  this  project  we  became  interested  in  the  possibility  of  developing  new  TL  materials  for 
temperature  sensing  using  the  SCS  method  (see  Section  3.2).  Attractive  points  of  this  technique 
include: 

®  Simplicity  and  cost-effectiveness; 

•  High  yield  (mass  content); 

°  Capability  to  achieve  high  purity  crystalline  materials; 

•  Capability  to  synthesize  multiphase  complex  oxide  powders. 

We  focused  on  demonstrating  that  the  SCS  method  has  the  potential  to  produce  materials  with 
high  TL  intensity  and  various  TL  peaks  that  could  be  useful  for  thermometry.  The  goal  is  to 
identify  potential  materials  and  understand  how  to  control  the  synthesis  conditions  to  engineer 
the  materials  with  the  desired  properties. 

Figure  54  shows  examples  of  X-ray  diffractograms  of  materials  synthesized  using  the  SCS 
technique,  demonstrating  the  possibility  of  obtaining  crystalline  samples. 
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Figure  54.  XRD  for  four  different  materials  produced  by  SCS:  La203,  A1203,  Ce02,  and  MgO. 


4.5. 1.1  Luminescence  from  SCS  samples 

As  discussed  in  Section  3.2,  several  compounds  were  synthesized  as  part  of  this  investigation. 
Figure  55  shows  examples  of  a  variety  of  ceramic  materials  prepared  by  SCS.  The  results  reveal 
a  variety  of  TL  curves  which  increase  with  the  dose  of  ionizing  radiation,  as  desirable  in  a  TL 
material.  In  these  examples,  the  TL  intensity  is  usually  not  very  high,  but  at  this  point  the 
experimental  parameters  for  synthesizing  these  materials  had  not  been  optimized. 
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Temperature  (°C)  Temperature  (°C) 

Figure  55.  TL  glow  curves  obtained  after  beta  irradiation  of  (a)  Al203:Er3+,  (b)  MgAl204:Ce3+,  (c) 
LaMgA103:Ce  ,  and  (d)  MgO.  All  samples  were  prepared  by  SCS  and  annealed  at  900  °C  during  2  h  in  air. 


An  important  part  in  engineering  new  materials  is  to  introduce  luminescence  centers  with 
emission  that  matches  the  detection  systems  used  for  TL  readout.  To  investigate  if  appropriate 
doping  can  be  used  with  the  SCS  technique  to  control  the  luminescence  centers,  we  synthesized 
MgO  doped  with  different  dopants  and  measured  the  RL  emission  spectrum  (Figure  56).  A 
tentative  assignment  of  the  emission  bands  in  the  RL  spectrum  observed  in  MgO  doped  with 
different  dopants  is  also  depicted.  The  emission  band  in  the  Ce3+-doped  MgO  is  attributed  to  the 
5d'4f°-»  4f'  transition  of  CeJ+  ion,  in  agreement  with  the  one  reported  by  Bapat  and  Sivaraman 
{,  1986  #683}  in  MgO:CeJ+,Li+.  Emission  from  Ce3+  is  known  to  vary  from  blue  to  the  red 
region  of  the  spectrum,  depending  on  the  crystal  field  (Blasse  and  Grabmaier,  1994).  The 
emission  from  Er  -doped  samples  show  lines  that  are  not  correlated  to  Er3+  emissions.  These 
emissions,  located  around  405  and  750  nm,  are  associated  with  Fe3+  and  Cr3+  impurities, 
respectively  (Takeuchi  et  al.,  1975).  In  Tm3  co-doped  sample,  the  sharp  emission  rays  around 
465  nm  are  assigned  to  'G4->  jH6  transition  of  Tm3+  as  it  has  been  reported  in  ZnS:Tm3+  and 
ZnO.Li  ,Tm  compounds  (Collins  et  al.,  1991;  Ronfard-Haret  and  Kossanyi,  1999),  while  the 
lines  around  360  nm  can  be  assigned  to  the  'D2^3H6  transitions  of  Tm3+  (Blasse  and 
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Grabmaier,  1994).  The  emission  lines  in  the  TbJ+-doped  samples  are  well-know  and  attributed  to 
the  indicated  transitions  of  the  Tb3+  ion  indicated  in  the  figure  (Bos  et  al.,  2006). 

These  results  indicates  that  that  Ce1+,  TmJ+,  and  Tb>+  are  entering  the  crystal  structure  and 
functioning  as  luminescence  centers  in  the  material. 


Figure  56.  RL  emission  spectrum  of  MgO  prepared  by  SCS  and  doped  with  different  lanthanides  (1  %  molar) 
and  Li  (10  %  molar).  The  samples  were  annealed  at  900  °C  during  2  h. 


Several  synthesis  parameters  have  been  investigated  to  control  and  optimize  luminescence 
processes  of  the  ceramics  produced  by  SCS,  including: 

°  Combustion  environment; 

°  Presence  of  co-dopants; 

°  Concentration  of  dopants. 

Figure  57  compares  A^O^iCe  and  AI2O3  samples  produced  by  SCS  carried  out  in  a  hot  plate 
and  inside  an  oven.  In  both  cases,  the  oven  combustion  resulted  in  higher  TL  intensity  than 
combustion  in  a  hot  plate.  We  infer  from  these  results  that  the  better  temperature  uniformity  in 
the  oven  increases  the  yield  of  ceramics  produced.  In  a  hot  plate,  the  temperature  gradients  are 


HDTRA1  -07-1  -0016 

University  of  Minnesota  (Twin  Cities) 


high  and  there  is  a  considerable  part  of  the  material  that  does  not  reach  a  temperature  sufficiently 
high  to  achieve  a  full  reaction. 

One  of  the  important  results  in  this  study  is  the  demonstration  that  Li+  co-dopant  can 
substantially  increase  the  TL  intensity  of  the  materials  being  synthesized.  As  an  example.  Figure 
58  shows  the  TL  curves  of  MgO:CeJ+  synthesized  with  or  without  the  presence  of  Li+.  The  Li+ 
co-doped  compound  shows  a  TL  intensity  which  is  1-2  orders  of  magnitude  more  intense  than 
the  sample  produced  in  absence  of  Li  .  In  addition,  new  TL  peaks  appear  in  the  temperature 
range  from  200-300°C.  It  is  worthwhile  to  note  that  MgO:Ce3+,Li+  materials  exhibit  several  TL 
peaks  in  a  large  interval  of  temperature,  as  it  is  intended  for  a  thermometric  application. 

Our  most  intense  TL  materials  prepared  so  far,  i.e.  MgO:Tm3+,Li+  and  MgO:Ce3+,Li+,  exhibit  a 
sensitivity  comparable  to  commercial  LiF:Mg,Ti  (TLD-100)  TL  dosimeters  (Figure  59). 
Although  the  main  TL  peak  of  Tnf+-doped  sample  is  located  at  too  low  temperature  for 
dosimetry  applications,  these  results  are  encouraging  since  they  demonstrate  the  potential  for 
producing  ceramic  materials  with  high  sensitivity  by  SCS. 

In  order  to  check  the  reproducibility  of  SCS,  two  samples  of  MgO:Ce+3:Li+  were  synthesized  in 
similar  conditions  and  their  TL  signal  was  then  acquired  (Figure  60).  Both  samples  show  the 
same  IL  peaks  with  only  a  slight  deviation  in  their  relative  intensity  and  their  temperature 
maximum  which  evidence  that  a  good  reproducibility  via  SCS  technique  should  be  achieved  in 
the  future. 

Finally,  quenching  curves  of  TL  and  RL  intensity  for  A1203  versus  Ce3+  concentration  in  the 
initial  solution  are  compared  in  Figure  61.  A  cerium  concentration  around  0.5%  is  the  optimum 
for  the  TL  intensity,  whereas  for  RL  the  optimum  concentration  is  around  0.2%.  This  difference 
between  the  two  quenching  optimal  concentration  is  related  to  the  fact  that  the  TL  process 
depends  not  only  on  the  concentration  of  luminescence  centers  (Ce3+),  but  also  on  the 
concentration  of  trapping  centers  in  the  material. 

The  possibility  of  controlling  the  trapping  centers  has  not  been  investigated  yet,  but  we  envision 
to  be  able  to  control  the  temperature  of  the  peaks  by  choice  of  suitable  co-dopants  introduced  in 
the  synthesis. 


HDTRA1-07-1  -0016 

University  of  Minnesota  (Twin  Cities) 


Figure  57.  Comparison  of  TL  from  Al203:Ce3+  and  Al203  powder  samples  produced  by  SCS  inside  an  oven 
and  in  a  hot  plate. 


Figure  58.  TL  glow  curves  obtained  after  beta  particle  irradiation  of  (a)  MgO:Ce3+(I%):Li+(i0  %)and  (b) 
MgO:Ce  ioj?%).  Both  samples  were  annealed  at  900  °C  during  2  h. 
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Figure  59.  TL  glow  curves  of  (a)  MgO:Tm3+(1%):Li+(I0  %>  annealed  powder  sample, 
sample  and  (c)  MgO:Ce3  (i%):Li+(10  %)  annealed  powder  sample  after  being  exposed  to 
(measured  with  10  mg  for  each  sample). 


Figure  60.  TL  glow  curves  obtained  after  beta  particle  irradiation  of  two  MgO:Ce: 
synthesized  in  the  same  experimental  conditions. 


(b)  TLD-100  powder 
0.2  Gy  beta  radiation 


(i%):Li+(io  %)  samples 


HDTRA1-07-1-0016 

University  of  Minnesota  (Twin  Cities) 


.8 

O 

T— 

£ 

I 

£ 

c 

S 

1? 

f 

C 


60 

50- 

40- 

30 

20-| 

10 

0 


(a) 


Batch  Number 


0 

#9 

A 

#10 

□ 

#26 

■ 

#27 

0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4 


Cerium  Concentration  (%) 


Cerium  Concentration  (%) 


Figure  61.  (a)  Integrated  TL  intensities  and  (b)  RL  intensities  of  AI203:Ce  samples.  Each  data  point  from 
batch  #27  represents  the  average  intensity  of  four  aliquots,  with  error  bars  displaying  the  standard  deviation 
of  the  mean.  Single  points  have  been  included  from  previous  combustions  (#9,  #10,  and  #26)  which  appear 
to  be  in  agreement  with  the  trend  set  by  the  behavior  of  batch  #27.  All  samples  were  given  a  dose  of  10  Gy 
prior  to  TL  measurements. 


4.6  Survivability  Testing  of  Luminescent  Particles  in  Explosions 

The  work  in  this  section  was  performed  in  collaboration  with  Jim  Lightstone  of  Indian 
Head  Naval  Surface  Warfare  Center. 

The  interior  of  an  explosion  is  one  of  the  harshest  sensing  environments  on  earth.  The 
violent  temperatures  and  pressures  will  destroy  all  conventional  sensors,  and  the  fireball  is 
usually  opaque  to  remote  optical  analysis  due  to  chemistry  or  debris  or  both.  However,  it  is  still 
critical  to  know  the  thermal  history  of  detonation  events  for  applications  such  as  the  sterilization 
of  weaponized  pathogens  and  military  explosive.  We  have  shown  earlier  in  this  report  that 
microparticles  of  thermoluminescent  oxides  may  be  used  to  measure  temperature  excursions  of 
hundreds  of  degrees  in  tens  to  hundreds  of  milliseconds.  Simulation  of  explosive  temperatures 
was  performed  using  microheaters,  but  these  lack  the  severe  pressure  and  shock  stimuli  which 
would  accompany  an  actual  explosion  and  so  therefore  could  not  assess  the  survival  of 
luminescent  microparticles  in  such  an  event. 

Thermoluminescent  (TL)  AI2O3.C  microparticle  sensors  are  demonstrated  to  survive 
detonation  of  pentaerythritol  tetranitrate  (PETN)  high  explosives.  Since  TL  particles  have 
previously  been  shown  to  act  as  sensors  for  rapid  high-temperature  events,  the  stability  of 
AbOsiC  particles  in  explosions  is  evaluated  and  several  conclusions  made.  First,  SEM  imaging 
reveals  particle  features  and  corners  with  radii  of  curvature  less  than  20  pm  yielding  to  rounded 
forms  and  increased  surface  roughness  post-explosion.  Chemical  analyses  by  XPS  and  confocal 
Raman  spectroscopy  show  small  increases  in  carbon  and  up  to  a  32%  decrease  in  aluminum  on 
the  particle  surfaces,  likely  indicators  of  A1  combustion.  Finally,  the  TL  emission  peak  and  shape 
remain  stable  to  within  1 .9%  and  4.0%  respectively,  less  than  experimental  error.  This  means 
that  the  particle  trap  density-of-states  and  luminescence  centers,  and  their  ability  to  preserve 
thermal  history  information,  are  largely  unaffected  by  detonation. 
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A  high  explosive  charge,  comprising  approximately  either  2  or  25  g  of  PETN  and  a 
detonator,  was  immersed  in  the  microparticles  and  then  triggered,  as  imaged  in  Figure  62. 
Particles  were  recovered  by  a  plastic  cone  placed  below  the  charge  assembly.  Chemical  analysis 
using  X-ray  photoelectron  spectroscopy  (XPS)  indicates  that  the  explosions  cause  a  relatively 
large  decrease  (Figure  63)  in  the  aluminum  present  on  the  particle  surfaces2;  coupled  with  a 
smaller  dip  in  oxygen  content,  this  suggests  the  combustion  of  atomic  aluminum.  A  small 
increase  in  carbon  is  also  observed.  Scanning  electron  microscope  (SEM)  examination  reveals 
that  the  AhC^C  microparticles,  originally  angular  bodies  with  sharp  comers  and  smooth  planes 
(Figure  64),  have  been  reshaped  into  much  more  rounded,  compact  configurations  with  high 
surface  roughness  (Figure  65).  Most  critically,  the  thermoluminescence  glow  curves  of  the  post¬ 
explosion  microparticles  were  collected  using  fresh  broad-spectrum  UV  irradiation  and 
compared  to  those  of  pristine  microparticles.  No  significant  difference  is  observed  among  the 
curves  (Figure  66);  the  full  width  at  half-maximum  and  peak  location  of  the  primary  A1203:C  TL 

peak  both  vary  by  less 
than  the  experimental 
error,  indicating  that 
the 

thermoluminescence 
characteristics  and 
mechanisms  of  the 
particles  have  not 
been  modified  by  the 
explosion  and  that 
their  ability  as  thermal 
history  sensors  has 
therefore  not  been 
Figure  62  -  The  detonation  of  the  charge  +  assembly  in  which  the  A1203:C  impaired, 
microparticles  were  embedded.  These  images  were  captured  as  20  ns 
exposures  at  5  ps  intervals. 


Binding  energy  (eV) 

Figure  63  -  High-resolution  XPS  spectra  of 
AfCbiC,  centered  on  the  primary  peaks  of  its 
elemental  constituents,  after  detonation  of  a  25  g 
PETN  charge. 


HDTRA1-07-1-0016 

University  of  Minnesota  (Twin  Cities) 


Figure  64  -  SEM  image  of  AI2O.VC 

microparticles  before  undergoing  high  explosive 
detonation. 


Figure  65  -  SEM  image  of  AfO^iC 

microparticles  after  detonation  with  25  g  of 
pentaerythritol  tetranitrate  (PETN). 


Temperature  (°C) 


Figure  66  -  Thermoluminescence  glow  curves  of 
AhCbiC  before  and  after  explosion  with  PETN 
charges.  The  data  were  collected  using  a  linear 
heating  ramp  of  approximately  0.9  °C/sec. 


4.7  Pulsed  Thermoluminescence 
4.7.1  Background 

Explosive  heating  creates  significant  differences  in  the  TL  compared  to  a  standard  slow 
heating  ramp.  This  section  will  examine  the  scientific  aspects  of  this.  There  are  a  number  of 
areas  of  research  that  could  benefit  from  the  capability  to  perform  TL  measurements  at  heating 
rates  orders  of  magnitude  higher  than  those  achievable  using  regular  TL  readers.  Examples  are: 

(i)  investigation  of  how  thermal  profiles  caused  by  detonations  affect  the  TL  of  microparticles; 

(ii)  development  of  techniques  to  perform  non-destructive  readout  of  the  TL  signal  for  radiation 
dosimetry.  This  is  similar  to  optically  stimulated  luminescence  (OSL)(Akselrod,  1999) 
(McKeever,  2003),  but  eliminates  the  need  for  optical  filters  used  to  discriminate  between 
stimulation  light  and  the  material’s  luminescence,  (iii)  basic  investigations  on  the  mechanisms 
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related  to  the  luminescence  efficiency  in  TL  materials;  and  (iv)  increasing  the  instantaneous 
intensity  of  thermo  luminescence  so  that  small  volumes  of  material  (e.g.  micro-  and  nano¬ 
particles)  will  emit  signals  above  detector  dark  noise  levels. 

Investigation  on  the  effect  of  temperature  profiles  generated  by  detonations  on  TL  of 
microparticles  could  lead  to  a  new  method  to  determine  the  thermal  profile  produced  by  these 
detonations  in  small  particles  or  organisms  such  as  bacteria  and  spores.  For  this  particular 
application,  contact  thermometry  and  stand-off  methods  are  not  able  to  determine  the 
temperature  profile  experienced  by  a  single  particle  during  the  entire  evolution  of  the  detonation 
event.  Luminescence  thermometry  has  been  demonstrated  using  nanoparticlcs,  but  a  practical 
application  has  yet  to  be  demonstrated  (Wang,  2002).  Since  the  TL  is  affected  by  the 
microparticle  temperature,  a  technique  to  determine  the  thermal  profile  could  be  developed  based 
on  the  comparison  of  the  TL  of  irradiated  microparticles  before  and  after  the  detonation. 
However,  this  requires  methods  to  experimentally  investigate  the  TL  behavior  at  the  fast  heating 
rates  experienced  by  the  microparticle.  For  example,  the  temperature  in  a  explosive-driven  shock 
wave  can  vary  by  thousands  of  degrees  Kelvin  in  a  few  microseconds  (Davis,  1 996).  Therefore, 
it  is  important  to  develop  capabilities  to  measure  the  TL  at  these  high  heating  rates  to  compare 
the  experimental  results  with  computational  models. 

Another  area  of  potential  interest  to  perform  fast  heating  rates  is  radiation  dosimetry.  Optically 
Stimulated  Luminescence  (OSL)  has  become  widely  adopted  in  dosimetry  due  to  several  factors. 
One  of  them  is  the  possibility  of  stimulating  the  crystal  in  a  very  controlled  way  by  using  fast 
pulses  of  light,  therefore  probing  the  OSL  signal  without  significantly  emptying  the  traps 
responsible  for  the  OSL  signal  (Akselrod,  1999).  Therefore,  fast  heating  and  cooling  rates  could 
also  lead  to  the  development  of  a  similar  approach  for  TL  dosimetry,  in  which  the  TL  signal  is 
probed  using  short  thermal  pulses  without  emptying  significantly  the  associated  traps,  as 
proposed  earlier.  One  of  the  advantages  would  be  to  eliminate  the  need  for  optical  filters  used  in 
OSL  to  discriminate  between  stimulation  light  and  the  material’s  luminescence. 

In  this  section,  we  describe  a  pulsed  thermoluminescence  (PTL)  technique  where  short  voltage 
pulses  are  applied  to  a  resistor  on  a  microheating  stage.  The  stage  heats  and  then  cools  rapidly, 
transferring  energy  to  one  or  more  luminescent  particles  on  the  stage,  exciting  trapped  carriers 
which  recombine  and  produce  light.  The  technique  is  demonstrated  using  particles  of  a-AhC^rC, 
which  is  a  common  TL  and  OSL  dosimetric  material. 


4. 7.2  Procedure 

Two  types  of  heaters  were  fabricated:  one  set  was  etch-released  from  the  substrate  and  had  slow 
time  constants  (~50ms);  the  other  set  was  not  etch-released  and  had  fast  time  constants  (less  than 
500ps).  These  latter  microheaters  consisted  of  silicon  nitride  plates  with  platinum  resistors  used 
for  joule  heating.  The  temperature  of  the  heaters  was  determined  by  measuring  the  change  in 
resistance  as  the  structure  was  heated. 

Platinum  was  chosen  for  the  resistors  since  the  temperature  coefficient  of  resistance  (TCR)  is 
constant  up  to  temperatures  near  600°C  (Firebaugh,  1998).  The  TCR  of  the  platinum  resistor  was 
determined  by  measuring  the  resistance  of  the  device  as  it  was  heated  on  a  hot  plate  up  to  90°C. 
The  resistance  versus  temperature  was  linear  over  the  entire  measurement  range  and  was 
assumed  to  be  linear  beyond  90°C.  We  calculated  the  temperature  change  as  a  function  of 
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applied  power  for  different  pulse  lengths  by  measuring  the  resistance  change  during  a  voltage 
pulse.  The  temperature  of  the  resistors  was  calculated  using  the  formula: 


(1) 


Where  AT  is  the  change  in  temperature,  R  the  resistance  during  heating,  Ro  the  resistance  at 
room  temperature  and  a  is  the  TCR  of  the  material. 

To  measure  the  luminescence,  the  photon  counts  were  recorded  with  a  photomultiplier  tube  for 
one  second.  This  time  period  was  chosen  since  it  was  much  longer  than  the  luminescent  lifetime 
of  Al203:C  and  thus  collected  essentially  all  the  photons  from  a  single  pulsed  excitation.  The 
temperature  pulse  occurred  immediately  after  the  photon  count  had  started,  which  allows  one  to 
assign  the  total  integrated  intensity  to  the  single  pulse. 

Immediately  after  this  measurement  was  recorded,  a  background  measurement  was  taken  for  one 
second.  This  verified  that  there  was  no  residual  luminescence  remaining  after  the  initial 
measurement.  Changes  in  this  measurement  were  never  observed  confirming  that  one  second 
was  long  enough  to  record  all  the  thermoluminescence  due  to  a  pulse. 

To  prepare  the  A^C^C  for  PTL  measurements,  a  single  particle  was  placed  on  a  microhcater. 
The  same  particle  was  used  for  all  measurements  and  was  approximately  80-  100pm  in  diameter. 
Its  size  was  estimated  by  comparing  it  to  the  size  of  the  heater  plate.  To  empty  the  traps  before 
exposure,  the  particle  was  heated  repeatedly  to  250°C  for  50  ms  until  the  measured  light 
intensity  was  equal  to  the  background  level.  The  traps  were  then  filled  by  exposing  the  particle  to 
205nm  light  for  20  min.  This  wavelength  coincides  with  the  absorption  band  of  F-centers  in 
AI2O3.  The  UV  photons  create  free  electrons  that  can  be  trapped  at  defects  in  the  crystal  (Lee, 
1979). 

There  were  two  different  methods  used  for  taking  PTL  measurements.  The  first  method  was 
temperature  ramping  PTL  (TR-PTL),  where  the  peak  temperature  of  each  pulse  was  increased  by 
a  constant  amount  over  that  of  the  previous  pulse.  This  allowed  the  measurement  to  be  displayed 
in  a  plot  of  photon  count  versus  temperature. 

The  second  method  was  constant  temperature  PTL  (CT-PTL).  In  this  technique,  the  temperature 
was  pulsed  repeatedly  at  the  same  level,  and  the  data  was  plotted  in  a  graph  of  photon  count 
verses  pulse  number.  This  allowed  the  change  in  trap  population  to  be  examined  as  a  function  of 
pulse  number. 

The  results  were  interpreted  according  to  the  first-order  TL  model,  in  which  the  probability  of  a 
charge  escaping  from  the  trap  per  unit  time  is  given  by: 


p  =  se 


(2) 


and  the  TL  intensity  I  is  given  by  the  rate  of  change  in  the  trapped  charge  population: 
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/  = 


(3) 


In  these  equations,  E  is  the  activation  energy  for  the  process,  s  is  the  frequency  factor  or  pre¬ 
exponential  factor,  k  is  the  Boltzmann  constant,  T  is  the  temperature,  and  n  is  the  concentration 
of  trapped  charges.  It  is  understood  that  first  order  kinetics  are  not  a  perfect  quantitative 
representation  for  AbOjiC  luminescence  because  of  quenching  at  higher  temperatures  (Akselrod, 
1999)  but  this  approximation  provides  a  simple  model  to  understand  the  qualitative  behavior  of 
the  system  and  demonstrate  the  potential  of  PTL.  The  intensity  as  a  function  of  temperature  for  a 
standard  TL  ramp  is  given  by: 


I(T)  =  n0se  h/kl  exp 


(4) 


where  p  is  the  heating  rate  and  T0  the  initial  temperature  of  the  sample.  The  TL  emitted 
during  the  /th  thermal  pulse  in  a  TR-PTL  or  CT-PTL  measurement  for  a  pulse  duration^  t  and 
initial  trap  population  of  n,  is  given  by  the  integral  of  (3): 

J  Idt  =  n  ( 1  -  exp  [- A/  •  se~,akTl  ] )  ^ 

The  surviving  concentration  of  trapped  charges  after  a  thermal  pulse  is  then: 
n,+]  =/7,exp[-At-i-e‘m7;]  ^ 

If 

1 »  exp  [-At  •  se~'Jk,‘  ] 

then  the  detrapping  is  small  between  pulses  and  the  intensity  remains  near  constant  for  all 
pulses. 


4.7.3  Data  from  Pulsed  Thermoluminescence  Experiments 

In  Figure  67  the  TR-PTL  curve  is  plotted  for  both  50  and  10  ms  pulses  as  curves  67b  and  67c 
respectively.  Each  point  corresponds  to  the  integrated  luminescence  measurement  from  one 
pulse.  Curve  67a  is  a  standard  TL  curve  of  a-AfeC^C  using  a  heating  rate  of  5  K/s  and  is  added 
for  comparison.  The  temperatures  listed  are  of  the  platinum  thermistor.  Note  that  the  TR-PTL 
peaks  have  a  similar  shape  to  a  standard  TL  curve  but  are  shifted  to  higher  temperatures.  The 
positions  of  the  peaks  were  seen  to  increase  with  decreasing  pulse  length.  Also,  the  temperature 
of  the  onset  of  luminescence  is  observed  to  increase  with  decreasing  pulse  length. 
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Figure  68  is  a  plot  of  CT-PTL.  The  curves  appear  to  decay  exponentially  with  pulse  number. 
However,  for  small  temperatures  the  decay  is  slow  enough  that  the  intensity  is  nearly  constant 
over  a  large  number  of  pulses.  Also,  the  intensity  was  observed  to  increase  with  the  temperature 
of  the  pulse. 


4.7.4  Discussion 

The  observed  temperature  shift  in  the  TR-PTL  curves  when  compared  to  the  standard  TL  curve 
has  two  causes.  First,  the  overall  amount  of  energy  delivered  to  the  heater  by  the  time  a  certain 
temperature  is  reached  differs  depending  on  the  heating  rate.  The  time  required  for  a  heater  to 
reach  a  temperature  of  500K.  from  room  temperature  at  a  heating  rate  of  1  K/s  is  twice  that  using 
a  heating  rate  of  2  K/s.  Clearly,  more  energy  is  delivered  to  a  heater  in  the  first  case.  This  means 
that  more  traps  will  be  excited  at  lower  temperatures  in  the  slow  ramp,  shifting  the  peak 
luminescence  lower  in  temperature.  This  dependence  of  the  TL  peak  position  on  the  heating  rate 
is  very  well-known  and  applies  equally  to  rapid  pulses.  A  TL  curve  obtained  using  a  heating  rate 
of  5  K/s  should  show  peak  luminescence  at  a  much  lower  temperature  than  a  TL  curve  produced 
using  a  series  of  10  ms  square  pulses  with  gradually  increasing  amplitude.  Likewise,  a  series  of 
similarly  increasing  50ms  square  pulses  should  lead  to  a  TL  curve  intermediate  between  the  two. 
This  is  the  behavior  measured  in  Figure  67  and  modeled  in  Figure  69  using  first  order  kinetics. 


Figure  67:  Curves  for  A1203:C  a)  intensity  of  a  standard  TL  curve.  Heating  rate  5°C/s.  b) 
Integrated  intensity  due  to  a  pulse  in  a  PTL  curve.  50  ms  pulses,  2°C  increase  per  pulse,  c) 
Integrated  intensity  in  a  PTL  curve.  10  ms  pulses,  2°C  increase  per  pulse. 


The  second  possible  reason  for  shifts  in  the  luminescence  is  due  to  the  imperfect  thermal  contact 
between  the  particles  and  the  heaters.  This  leads  to  a  delay  in  heat  transfer  from  heater  to 
particle.  If  the  delay  is  longer  than  the  pulse  time,  the  particle  would  never  able  to  reach  the 
expected  temperature.  One  may  be  seeing  the  beginnings  of  the  effect  by  the  larger  separation 
between  the  50  ms  and  10  ms  peaks  in  the  measured  data  over  the  simulations.  Based  on  the 
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close  match  between  Figures  67  and  69,  this  appears  to  be  a  secondary  effect  that  could  be 
eliminated  by  using  deposited  thin  films  on  microheaters  as  luminescent  materials. 

The  potential  use  of  PTL  to  allow  a  dosimetric  material  to  be  read  repeatedly  can  be  seen 
in  Figure  68.  At  lower  temperatures  CT-PTL  is  near  constant  over  many  pulses.  This  indicates 
the  luminescence  can  be  read  multiple  times  from  the  same  particle  without  losing  stored 
information  about  the  population  of  traps. 


100  150  200  250  300 

Pulse  Number 


Figure  68:  PTL  at  constant  temperature.  The  pulse  length  was  50  ms.  The  temperatures  of  the 
pulses  for  different  curves  were  175, 180, 185, 190, 195  and  200°C. 


Figure  69:  Simulation  first  order  curves  for:  a)  TL  5  °C/s  b)  TR-PTL  50ms  c)  TR-PTL  10ms.  This 
data  is  not  expected  to  exactly  match  that  of  Figure  2  because  a-ALC^C  is  not  an  ideal  material 
with  first-order  kinetics,  but  the  qualitative  shift  is  seen  to  be  very  similar. 
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The  use  of  PTL  may  also  be  a  helpful  technique  in  decreasing  the  effects  of  thermal  quenching  in 
luminescent  particles  as  well  as  potentially  investigate  the  thermal  quenching  mechanism.  For 
this  to  be  a  viable  option,  the  thermal  time  constant  of  the  microheater-plus-TL-material  would 
have  to  be  much  less  than  the  F-center  lifetime  (35  ms  at  room  temperature)  (Akselrod,  1999) 
Then  by  applying  temperature  pulses  shorter  than  the  F-center  lifetime,  the  majority  of  the  decay 
can  occur  at  room  temperature  where  the  signal  is  not  affected  by  thermal  quenching.  Since  the 
heaters  in  this  study  have  time  constants  on  the  order  of  500ps,  these  studies  will  be  possible  by 
matching  the  thermal  mass  of  the  TL  material  with  improved  thermal  contact,  such  as  would  be 
obtained  using  deposited  thin  films. 

The  technique  of  PTL  could  also  be  used  in  an  integrated  badge  dosimeter  which  would  be  able 
to  warn  the  user  if  they  had  received  a  harmful  amount  of  radiation.  The  device  could  consist  of 
a  TL  particle,  microheater  and  avalanche  photodiode.  Since  an  avalanche  photodiode  tends  to  be 
noisier  than  a  PMT,  a  single  PTL  pulse  could  potentially  be  used  to  increase  the  intensity  of 
emitted  light  above  the  noise  level  of  the  photodiode  in  a  single  temperature  pulse.  Although  in 
our  work  we  were  unable  to  test  this  due  to  equipment  limitations,  this  can  be  anticipated  to 
increase  the  intensity  of  the  signal  by  imagining  a  pulse  at  the  same  temperature  as  the  peak  of  a 
standard  TL  curve.  The  luminescent  efficiency  and  rate  of  trap  emptying  would  be  the  same  in 
both  cases.  However,  with  a  single  temperature  pulse,  there  was  no  previous  emptying  of  traps  at 
lower  temperatures  like  there  would  be  with  a  standard  TL  ramp.  Since  there  would  be  a  greater 
number  of  filled  traps,  the  intensity  is  anticipated  to  be  greater. 


5  CONCLUSIONS 

Construction  of  a  thermoluminescence  system  with  controlled  explosive-rate  heating  and 
cooling: 

•  Microheaters  for  simulating  explosive  temperatures  in  microparticles  have  been 
completed.  Various  sorts  of  microhcaters  have  fabricated.  The  fastest,  for  rapid  heating 
can  ramp  to  350°C  with  time  constants  of  a  few  hundred  microseconds  with  equally  fast 
cooling  rates.  Other  heaters  have  been  fabricated  with  maximum  temperatures  that  can 
exceed  800°C  but  with  time  constants  in  the  tens  of  milliseconds. 

®  A  thermoluminescence  system  has  been  constructed,  tested,  and  verified  that  simulates 
explosive  heating  without  confounding  pressure  effects.  The  entire  process  from  UV 
excitation  for  trap  filling  of  microparticles  to  data  processing  of  TL  data  has  been  fully 
automated  for  both  standard  TL  and  explosive  heating. 


Demonstration  of  basic  thermometry  concept: 


o  Basic  modeling  to  extract  both  maximum  temperature  and  cooling  time  has  been 
performed.  Modeling  of  the  luminescence  with  ramp  speed  and  pulse  energy  was 
completed.  Simulations  show  that  obtaining  the  peak  temperature  that  a  particle  has 
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reached  is  relatively  straightforward  in  theory  but  information  on  cool-down  is  more 
difficult  but  not  impossible.  The  important  parameters  of  the  microparticles  are  trap  depth 
and  lattice  (phonon)  interaction  time.  Using  multiple  particles  is  necessary  in  principle  to 
isolate  maximum  temperature  and  cooling  time. 

o  Using  experimental  and  numerical  studies  based  on  LiF:Mg,Ti,  it  was  demonstrated  that 
a  superposition  of  first-order  TL  peaks  offers  more  information  than  a  single  TL  peak. 

We  showed  that  two  different  profiles  that  would  result  in  the  same  change  of  a  single  TL 
peak  (e.g.  a  fast  heating  to  a  high  temperature  versus  a  slow  heating  to  a  lower 
temperature)  would  affect  differently  other  TL  peaks  in  the  material,  therefore  opening 
the  possibility  of  distinguishing  the  heating  profiles  based  on  the  superposition  of  TL 
peaks. 

®  The  core  concept  of  the  program  of  using  microparticle  TL  to  extract  maximum 
temperature  with  rapid  heating  rates  and  slower  cooling  times,  mimicking  what  would  be 
seen  in  an  agent  defeat  test.  MgSi04:Co,Tb  and  ALOsiC  microparticles  were  used  with 
microheater  excitation.  Peak  temperatures  were  recovered  with  good  accuracy.  The  peak 
temperatures  ranged  from  232  °C  to  313  °C  with  durations  less  than  200ms.  The 
observed  change  in  the  thermoluminescence  peak  intensity  ratio  of  two  Mg2Si04:Tb,Co 
traps  has  been  predicted  with  an  average  error  of  4.4%  using  first-order  kinetics 
modeling. 

o  In  the  process  of  creating  explosive  heating,  we  discovered  a  new  characterization 
technique,  pulsed  thcrmoluminescence,  that  that  enables  one  to  assess  the  trap  population 
using  a  single  or  series  of  pulses.  This  may  allow  single  particles  to  have  an  instantaneous 
thermoluminescent  intensity  that  is  higher  than  standard  TL,  enabling  the  read-out  of 
very  small  particles,  which  is  particularly  important  for  luminescent  thermometry. 


Characterization  of  existing  dosimetric  materials  for  thermometry: 

®  Existing  materials  are  available  with  properties  that  could  be  used  as  a  temperature 
sensor.  LiF:Mg,Ti  emerged  as  the  most  favorable  material  due  to  its  high  TL  intensity, 
presence  of  several  TL  peaks,  and  agreement  between  experimental  data  obtained  using 
this  material  and  simulations  using  TL  models. 

©  One  of  the  predominant  problem  for  the  use  of  existing  TL  materials  in  temperature 
sensing  is  the  light  sensitivity  of  the  trapping  centers.  Except  for  LiF:Mg,Ti  and  natural 
topaz,  all  other  materials  show  some  degree  of  light  sensitivity. 

o  Investigation  of  new  nanophosphors  produced  by  Solution  Combustion  Synthesis  (SCS) 
indicate  a  possible  new  route  for  material  synthesis.  The  oxyorthosilicates  investigated  in 
this  study  were  not  suitable  for  temperature  sensing  because  of  the  low  TL  intensity,  but 
this  does  not  exclude  the  possibility  of  using  SCS  to  produce  other  materials  with 
desirable  properties. 
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Development  of  new  materials  for  luminescence  thermometry: 

o  The  possibility  of  synthesizing  TL  materials  using  the  Solution  Combustion  Synthesis 
(SCS)  technique  was  demonstrated.  Materials  were  produced  with  favorable  TL 
properties,  i.e.,  several  TL  peaks  and  emission  bands  corresponding  to  the  emission  from 
dopants  introduced.  However,  this  is  a  new  field  of  development  in  TL  research  which 
will  require  substantial  effort  before  the  materials  can  be  properly  engineered  with  the 
desired  properties. 
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